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CHAPTER 1 


DISTURBING TORQUES AND POINTING ERRORS 
FOR LARGE EARTH-ORIENTED MICROWAVE REFLECTORS 


by 

Richard H. MacNeal 




PREFACE 


The successful performance of any structure depends largely 
on the identification of the critical or primary loads eind design 
criteria on which the design is based. The future promises large 
structures which must be deployed, erected, assembled, or fabri- 
cated in space. For such structures, which will not be required 
to face the launch environment, the primary design requirements 
will be derived from the space-flight environment and will deal 
with phenomena as primary criteria which have been considered as 
only secondary in the past. The design of such genuine "space" 
structures will require a solid foundation of critical criteria. 


Astro Research Corporation has a contract from NASA Langley 
Research Center to study critical design criteria for large space 
structures. The objective is to identify and establish critical 
baseline design requirements for a family of structures by a series 
of rational parametric analyses. The results, presented clearly 
and in detail, should improve the basis for future space structures 
system and technology efforts. They will also form the beginning 
of the needed solid foundation of design criteria. 

This report is one of a series dealing with critical design 
requirements for large space structures. In particular, this 
chapter investigates the disturbing torques on large Earth-oriented 
P^J^3,bolic reflectors and of methods to control the resulting point- 
ing errors. The chapter considers an altitude range and ranges 
of lector size, mass density, focal ratio, and pointing error 
which are sufficiently broad to cover most applications. This 
chapter was prepared at MacNeal— Schwendler Corporation under a 
subcontract from Astro Research Corporation. 
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LIST OF SYMBOLS 


Root Symbols 
A = area 

B * magnetic field intensity 
C = dimensionless coefficient 
D = diameter of reflector 
f = focal distance 
g = damping coefficient 

G = acceleration due to gravity at Earth's surface 
h = altitude 
H = angular momentum 
I ■ moment of inertia; current 
Igp “ specific impulse 

k = dimensionless 'coefficient 
S. = distance 

m » mass 

m' ” mass per unit area of reflector aperture 
M =» moment 
p = pressure; d/dt 
r - radius 

R » distance to center of Earth 
R' = resistivity (ohm-meter) 
t = time 

T = tension 

V = velocity; volume 

X = distance from c.g. in forward direction along flight path 


y = distance to right from c.g., looking forward 
z = distance from c.g. downward toward Earth 
a = angle of attack; dimensionless coefficient 

Y = angle between sunline and normal to surface (or pointing axis) 

e = small number 

p » density 

o “ stress 

(j> = roll angle 

9 = pitch angle 

ij; = yaw angle 

u> = frequency, rad/ sec 

ft = orbital angular velocity, rad/sec 

Subscripts 
a , anchor 
c , control 
d , drag 

e , expelled mass 

0 , steady value; overturning moment 

1 , inertia 

m , moment; mission 
r , reflector 

s , solar radiation pressure 
w , inertia wheel 


SECTION 1 
INTRODUCTION 

One of the earliest types of large space structures will be parabolic 
microwave reflectors with diameters ranging from 100 meters to 1000 meters. 

An important application of such reflectors will be to conduct surveys of 
the Earth’s surface in order to obtain prompt and detailed information re- 
garding oceanographic and biological phenomena. 

This chapter is part of a study of large space structures undertaken 
by Astro Research Corporation for NASA under Contract NASl-15347. It reports 
the results of a general investigation of the disturbing torques on large 
Earth-oriented parabolic reflectors and of methods to control the resulting 
pointing errors. The chapter considers an altitude range and ranges of 
reflector size, mass density, focal ratio, and pointing error which are 
sufficiently broad to cover most applications. 
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SECTION 2 

MISSION PARAMETERS AND CONTROL SYSTEM FUNCTION 

The range of mission parameters to be investigated is shown in Table 1. 

The function of the control system is to maintain the pointing accuracy within 
the range of tolerance indicated in Table 1 (10~^ to 10”^ radian). No tolerance 
is given for the angular deviation about the pointing axis (yaw direction), but 
it can be assumed that some tolerance must exist, if the Incoming radiation is 
polarized, or if the microwave beam is to be electronically scanned. 

Since the orbit is assumed to be circular, the nominal pointing direction 
is toward the center of the Earth (geotropic). Thus, the vehicle must rotate 
at constant angular velocity once per orbit in the pitch direction. No excur- 
sions of any kind from this pointing direction are desired. 

It will be noted from Table 1 that the ranges of the mission parameters 
and the pointing error tolerances are very broad. As a result, no single 
type of control system is optimum over the entire spectrum of p>arameters 
and tolerances. Unconventional control systems should be considered because 
the size and mass density parameters are outside the current state-of-the-art. 
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SECTION 3 

OVERTURNING MOMENTS AND POINT ERRORS DUE TO ENVIRONMENTAL EFFECTS 

The principal potential sources of pointing error in the vehicle's 
operational environment are 

1. Temperature Gradient and Thermal Shock 

2. Atmospheric Drag 

3. Solar Radiation Pressure 

4. Micrometeoroids 

5. Gravity 

Each of these items is discussed in separate subsections. Figure 1 
summarizes and compares equivalent loading pressures as a function of orbital 
altitude due to atmospheric drag, solar radiation pressure, and the combined 
effect of gravity gradient and unsymmetrical mass distribution. The follow- 
ing conditions are assumed in Figure 1. 

• Drag Coefficient, C^ = 2.5 

• Normal incidence and perfect reflection of solar radiation 

•Worst case orientation (45“) of the mass unbalance (mass per 
unit area times distance) from the geocentric direction, i.e., 
the pointing axis. 

Within the range of orbital altitudes considered in this study (200 km to 
geosynchronous), one or another of the three effects dominates the others 
at different altitudes. 

Maneuvers are not a source of pointing error, because it is not re- 
quired that the vehicle exercise maneuvers. Moments due to interaction 
with the Earth's magnetic field, which are significant for some spinning 
spacecraft, are not considered to be significant for the mission under 
consideration. They can, however, be used as a means of active control, 
see Section 4,6. 
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3. 1 Temperature Gradient and Thermal Shock 


Temperature gradients occur for various reasons including the facts 
that some parts of the vehicle are shaded from the Sun by others, and that 
the entire vehicle may pass into the shadow of the Earth. The main effect 
of temperature gradient on pointing accuracy is to cause thermal strains 
which distort the shape of the reflector and thereby change the pointing 
axis. Such errors cannot be detected by simple devices, and they are 
beyond the scope of the present investigation. Appendix A includes a very 
simple analysis of the effect on pointing accuracy of thermally Induced 
length changes of the feed supports. It is shown that the pointing errors 
in radians is of the same order of magnitude as the differences in the 
thermal strains, and that the allowable differential strain decreases as 
f/D increases. 

The effects of thermal shock caused by a rapid change in temperature 
are similar in nature, except that inertia forces are brought into play. 
Thermal shock occurs when the vehicle passes into the Earth's shadow. At 
500 km altitude, for example, the time to cross the penumbra is about 8 
seconds, so that structural modes with periods this long or longer will 
respond dynamically. Cumulative changes in the pointing direction due to 
nonlinear kinematic effects will be much smaller than the linear effects 
and may be disregarded, provided that some form of directional stability 
is provided. 
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3.2 Aerodynamic Drag 


Newtonian flow theory may be used to estimate the aerodynamic forces 
and moments at orbital altitudes. In Newtonian flow, the pressure is normal 
to the surface and has the magnitude 

P = • sin^a (1) 

where 

1 2 

-^pV = dynamic pressure 

Cjg *= normal force coefficient (2. 5) 
a = angle of attack 

Both density and velocity are fvinctions of altitude. The following table 
gives values of p vs. altitude for a = tt/2. 


h 

(km) 


300 

400 

500 

600 

700 

800 

900 

1000 


h - y 

N/m^ 

300x10"^ 

60 

20 

6 

2 

.8 

.4 

0.15 


It is noteworthy that is less than Pg, the solar radiation 
pressure at normal incidence, for altitudes greater than 800 km, see 
Figure 1. 

If the vehicle is symmetrical, the only component of moment due to 


aerodynamic pressure is in the pitch direction (normal to the plane of the 
orbit). Incidental or accidental lack of symmetry causes small moments in 
the roll and yaw directions. 
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If the surface of the reflector is smooth and if it is assumed, for 
convenience, to be spherical, then, clearly, the center of pressure passes 
through the center of the sphere. The moment arm between this point and 
the center of gravity is slightly less than 2f, where f is the focal distance. 
An expression for the moment is 

“d • ?d ( 2 ) 

where D is the diameter of the reflector and is a coefficient which 
depends on the focal ratio as shown in the following table 

m 

0.5 .0117 

1.0 .00277 

2.0 .000672 

The values of are derived in Appendix B, If the surface of the reflector 
is not smooth, (for example, a surface consisting of a grid work of wires) 
the force per unit area will be greater, but the moment arm will be less, 
resulting in a moment coefficient of the same order as in the table shown 
above. It is also possible, with non-smooth surfaces, to balance the drag 
forces by design so as to achieve zero net pitching moment. 

In order to appreciate the magnitude range of the pitching moment due 
to aerodynamic drag, consider the following table which records vs. D and 
f/D for h = 500 km. 



100 m 300 m 1000 m 
M (Newton-meters) 


0.5 

1.9 

49.7 

1842 

1.0 

0.5 

11.7 

434.5 

2.0 

0.1 

2.9 

106.1 
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For h » 300 km, the moments are an order of magnitude larger and for 
h » 700 km, they are an order of magnitude smaller. For non-solid surfaces, 
the moments should be reduced by the surface solidity ratio. 

Since the orientation of the vehicle remains constant with respect to 
the flight path, the dynamic component of aerodynamic moment is very small. 
Control surfaces may, however, be used to produce dynamically varying aero- 
dynamic moments, see Section A. 4. 
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3.3 Solar Radiation Pressure 


Solar radiation pressure is a significant source of overturning moment 
for large spacecraft. For a perfect reflector, the pressure is normal to 
the surface and has the magnitude 

p = Pg cos^Y (3) 

where Pg = .9x10”^ N/m^ near the Earth 

Y “ angle between Sun line and normal to surface. 

For a perfectly absorbing surface which radiates equally from both sides, 
the pressure is directed away from the Sun and has the magnitude 

p - -i pg COSY (4) 

Equation (3) will be used for the purpose of estimation, except that pg 
will be reduced to ,75x10”^ N/m^ to account for absorptivity. 

A method has been developed in Reference 1 to calculate the moment 
due to radiation pressure on a perfectly reflecting cylindrical segment 
of small Included angle. The moment per unit length about the midchord 
of the strip, as predicted by the method, is 

“ 12 sin2Y/Rc (5) 

where C is chord of the cylindrical segment, is its radius of curvature, 
and Y is the angle between the normal to the surface and a line pointing 
at the Sun* Note that the magnitude of the moment is maximum for y ® ^5® 
plus multiples of 90^* For Sun angles near 90® and 270®, some parts of the 
reflector shade other parts and Equation (5) becomes inapplicable, although 
it is clear, in any case, that the moment should change sign near y = 90®* 
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Equation (5) is adequate to make rough estimates of the moments on a 


spherical reflector, with the following identifications: 

C = D, the diameter of the spherical reflector 
= 2f , the radius of the spherical reflector 

c IT p2 , the moment per unit area 


then 

where 

''s " 2|-f 

The following table evaluates for various values of the focal ratio, 
f /D, and y = ^5° . 

f/D tcp 

0.5 .0833 

1.0 .0417 

2.0 .0208 

Since the axis of the reflector points toward the center of the Earth 
at all times, the angle between the Sun line and the axis will vary in a 
manner that depends on the orbit. If the Sun line is normal to the plane 
of the orbit, then y = 90® at all times, and the overturning moments are 
small and steady. If the Sun line is tangent to the plane of the orbit, y 
will rotate through 360® in one orbital revolution and the pitching moment 
will exhibit a strong 2/rev component. Such trajectories may also pass 
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through the Earth’s shadow, further complicating the harmonic content of 
the overturning moment. The transfer from light to darkness through the 
penumbra at an altitude of 500 km requires about eight seconds which is less 
than one-thousandth of the orbital period. The gradient of light intensity 
across the reflector disk is not significant during passage of the penumbra, 
which has a length of 60 km at 500 km altitude. In addition, the angle be- 
tween the Sun line and the plane of the orbit will show an annual variation, 
unless the orbit is in the ecliptic plane. 

For non-solid smooth surfaces, the moments should be reduced by the 
surface solidity ratio. For a surface consisting of wires with circular 
cross-section, the moments will be zero unless some parts shade others. 

For real structures, the calculation of k in Equation (6) will require 

s 

detailed analysis. 
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3.4 Micrometeoroids 


A worst case analysis of the pointing error due to meteoroid impact 
is worked out in Appendix C. It is shown there that the maximxmi error 
angle due to multiple impacts of small meteoroids (micrometeoroids) is 


®max 


k(m*)^ 

D2(m3g)l/2 


( 8 ) 


where k is a constant, m is the natural frequency in pitch, and g = 2C/Cj, 
is the damping coefficient of rigid body pitching oscillations. 

In Eqtiation (8), it is assumed that the reflector is a shell of uni- 
form thickness which is penetrated by large micrometeoroids. Open con- 
struction,. or concentrated reinforcement, increases the probability of 
severe structural damage and occasional large pointing disturbances due 
to heavy impacts. 

If the following worst case values are assumed: 

®max “ radian 

m' “ 1 kg/m^ 

D = 100 m 

u) “ 10 rad/sec (approximate orbital frequency 
for altitudes below 1000 km) 


then the required amount of damping is g = 7.3x10“^. This is a very small 
value which should be exceeded for other reasons. 

Although micrometeoroids do not appear to be a significant source of 
pointing error for the range of vehicle parameters considered in the present 
study, it is clear from Equation (8) that their effect is magnified for the 
smaller, denser vehicles which typify conventional spacecraft design. 
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Also, the possibility of severe structural damage due to meteoroids 
should not be overlooked, particularly for vehicles with a large running 
length of structural reinforcement. Redundant design is the order of the 
day for such members (see Reference 2). 
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3.5 Gravity Gradient 


Gravity gradient is a strong source of overturning moment for orbiting 
spacecraft which are required to point in arbitrary directions. In the 
present case, the desired pointing axis is toward the center of the Earth, 
and the vehicle is in a circular (or nearly circular) orbit. For this case, 
the vehicle can be designed so that gravity gradient is a source of stability 
rather than a source of disturbing torque. If the principle axis of inertia 
with the least polar moment of inertia is aligned with the pointing axis, 
then gravity gradient will tend to resist overturning moments (see Appendix 
D). If, in addition, the intermediate principle axis of inertia is aligned 
with the flight path, gravity gradients will tend to resist yawing moments 
(moments about the pointing axis) . 

If the principle axes of inertia are not aligned with the pointing axis 
and with the flight path, gravity gradient will produce moments whose magni- 
tudes may be computed by Equation (6) of Appendix D. They are proportional 
to the cross-products of inertia which are under design control (within manu- 
facturing tolerance). 

A more serious requirement is that the pointing axis be the axis of 
least polar moment of inertia in order to achieve stability. This is 
illustrated by the following table which indicates the minimum ratio of 
feed mass to reflector mass required to obtain stability for reflectors 
with various f/D ratios. (The mass of the reflector is assumed to be uni- 
formly distributed and the feed supports are assumed to be massless.) 

°feed/”S:efl. 

0.5 .50 

1.0 .07 

2.0 .015 
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Since the indicated feed mass may be larger than that desired for 
other reasons, a potential stability problem arises which can be solved 
either by placing additional mass on the pointing axis or by providing an 
active control system. This subject is treated later. 

Another important effect is that, if the orbit is noncircular, the 
pitching angular velocity required to point the reflector toward the center 
of the Earth is not constant. This effect, although Important, is beyond 
the scope of the present investigation. It results in a scanning velocity, 
as measured at the surface of the Earth, that is nonuniform, but predictable. 

The effects of variations in the magnitude of the gravitational constant 
(oblateness) will be discussed later in connection with the stability of a 
gravity-gradient control system. Section 5.1. 
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SECTION 4 

CANDIDATE CONTROL SYSTEMS 

As applied to attitude control, a control system consists of a means 
to detect angular error, a means to apply a mechaniccal moment to correct 
the error, and a means to modify the moment in response to the error. A 
control system is passive if all three elements of the system are combined 
in the same mechanism. Otherwise, it is active. The advantages usually 
attributed to passive control systems are that they are more reliable and 
less subject to instability. Their principle disadvantage is larger weight 
for the same error . 

Emphasis will be placed, in the present discussion, on the means to 
apply a moment. The other elements of an active control system are con- 
ventional and have no unusual requirements in the present application. 

The source of moment to correct angular error may either be internal 
(expelled mass or stored angular momentum^ or it may be external, deriving 
from environmental effects such as gravity gradient and solar radiation 
pressure. The external sources of moment are particularly interesting for 
the present application because their intensity is weak compared to internal 
sources, and consequently they require large physical dimensions in order to 
be effective. For the same reason they are not particularly attractive for 
small spacecraft. 

The properties of seven different sources of control moment are 
summarized in Table 2. The first three employ internal sources of control 
moment and the last four employ external (environmental) sources. Formulas 
for the control moments provided by each of the seven sources are given in 
Table 3. Each of the candidate control systems will be discussed in separate 
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subsections, but it is useful to examine the summary information presented in 
Tables 2 and 3 before going into the details. 

It will be noted from Tables 2 and 3 that the internal sources of 
control moment are not altitude dependent, in contrast to the external 
sources which are, with the exception of solar radiation pressure which 
depends only on distance from the Sun, The effect of radiation pressture 
is, however, so weak that it is masked by aerodynamic pressure at low 
altitudes, as Indicated in Table 2, The same is true of gravity gradient 
and the Earth's magnetic field with the added restriction of a practical 
upper altitude limit due to 1/R^ decay of both effects (see Table 3, last 
column), Aerod3mamic pressure decays rapidly with altitude and is unusable 
as a source of control moment above 1000 km altitude. Thus, it is seen 
that no single external source of control moment is practical within the 
entire altitude range of the present study (200 - 30,000 km). 

The practical size ranges shown in Table 2 are derived from detailed 
analysis supporting the present study, but they can be inferred, to some 
extent, from the formulas presented in Table 3 and from previous discussion. 
The discussion of Section 3 showed that aerod3mamic pressvire and solar 
radiation pressure are the largest sources of overturning moments for the 
present application. Therefore, the use of the same sources to provide 
control moments should not result in any preferences with regard to size. 

The control moment available from gravity gradient, on the other hand, is 
seen to depend on the fourth power of vehicle size rather than on the third 
power. Thus, gravity gradient should be a powerful source of control moment 
for large vehicles, as indicated in Table 2, 
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The instantaneous value of the control moment available from expelled 
mass can be very high, but the average long term value is limited by weight 
considerations. Thus, the major limitation on the use of expelled mass (in 
particular, chemical propulsion) is the duration of the mission, which must 
be long for large expensive spacecraft in order to return initial invest- 
ment# The size of the vehicle does not directly affect the effectiveness 
of expelled mass relative to aerodynamic pressure and solar radiation 
pressure. 

Control systems which utilize stored angular momentum are penalized 
in the present application by a limitation on their size and by the constant 
vehicle angular velocity required by geotropic pointing. The size limita- 
tion derives from the fixed diameter of the launch vehicle and from the 
assumption that it will be impractical, for the foreseeable future, to 
construct a high performance inertia wheel in space. Fixed size degrades 
performance relative to environmental effects because the stored angular 
momentum, H, is proportional to the diameter of the wheel , given a fixed 
mass and optimized material selection, whereas the overturning moments 
are proportional to the cube of the diameter of the vehicle . 

Geotropic pointing limits the axis of the wheel to be the pitch axis 
(normal to the orbital plane) with the result that control from stored 
angular momentum can only be achieved by a combination of momentum bleeding 
(pitch axis) and gyrostabilization (roll and yaw axes). Gyrostabilization 
is particularly weak for large vehicles as shown by the formula in Table 3. 

The conclusions in Table 2 regarding use of the Earth’s magnetic field 
are based on a study. Reference 3, in which the vehicle was very large 
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(1500 m) and very light (.001 kg/m^). For this case magnetic field inter- 
action appeared to be a very effective means of control. 

The answers to the questions regarding passive and active control 
posed in Table 2 are in most cases self-evident, particularly if the re- 
quirements of the present application are considered. For example, a 
passive control system based on aerodynamic pressure (i,e,, simple aero- 
dynamic control surfaces) is possible because the velocity vector has a 
fixed orientation with respect to the vehicle. Solar radiation pressure 
cannot be so used because its direction and magnitude change continuously 
relative to vehicle axes. 
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4,1 Expelled Mass 


If gas jets are placed at four equally spaced points on the perimeter 
of the reflector, then the average available long-term restoring moment 
about a transverse (pitch or roll) axis is 


where 


G » acceleration of gravity 
m^ = stored mass of expellant 
I-_ = specific impulse 

b P 

D = diameter of reflector 
= mission duration 


G m I 

— 


m 


D 


(9) 


In order to estimate the required value of m^, assume that opposes 
an overturning moment due to solar radiation pressure. Solar radiation 
pressure provides a good example because it is independent of orbital 
altitude. Equations (6) and (7) of Section 3.3 provide an expression for 
the maximum overturning moment due to radiation pressure. If it is 
assumed that f/D = 1.0, and that the average overturning moment is equal 
to one-half of the maximum value, then the overturning moment is 


% = .0208 ‘^PgD^ 

In Equation (9), substitute 


( 10 ) 




m' 

e 


• - 2 - 1)2 


( 11 ) 
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Then, by equating to 

( 12 ) 

which shows that, m’ , the propellant mass per unit of reflector area, is 

e 

independent of vehicle size. As an example, select 

Pg » .75x10”^ N/m^ 
tjjj = 10 years = 3.15x10® sec 
G = 9.81 m/sec^ 

Igp = 200 sec 

where the value of I^p is typical for chemical propulsion. 

For this example, = .1002 kg/m^. Referring to Table 1, the range 
of the mass per unit area of the reflector is from .01 to 1.0 kg/m^. It 
is seen that the required propellant mass is excessively large, at least 
for the low range of reflector mass density. 

The required propellant mass to overcome solar radiation pressure can 
be reduced dramatically by using open wire construction for the reflector 
rather than solid construction, as mentioned in Section 3.3. 

Much lower values of m' can be achieved by using solar electric pro- 
pulsion (ion propulsion) which electrically accelerates charged particles. 
Energy is provided by solar cells. Values of I^p of the order of 10,000 
are technically possible for ten years of continuous thrust. Reference 4, 
but the level of technical complexity is high. 
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At low altitudes an exactly similar method may be used to estimate 
the expelled mass required to resist overturning moment due to aerodynamic 
drag. A related problem, which may set a minimum feasible altitude, is 
orbital decay due to aerod3mamic drag. Mass can be expelled continuously 
at a low rate to resist the drag force* The force exerted by expelled mass 
is 


e sp 


m 


(13) 


The drag force for a smooth reflector is 


The drag coefficient is evaluated in Appendix B as a function of f/D. 
For f/D “1.0, “ .00138. Using this value, the required expelled mass 

per unit area is, from Eqxaations (11), (13) and (14) 

m' “ _:.00138 p^ 

® ^ ^sp 

As an example, select 

m’ =* .01 Kg/m^ 

8 

tjjj = 10 years = 3.15x10® sec 
G = 9.81 m/sec^ 

Igp - 200 sec 

and calculate and the corresponding altitude. For this example p^ * 
4.51x10 ^ and the altitude is, from the short table given in Section 3.2, 
about 650 km. Increasing the expellant mass by a factor of ten decreases 
the altitude to about 450 km. The use of solar electric propulsion would 
appear to be the only practical method to achieve even lower altitudes for 
the application considered in this report. 
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4.2 Stored Angular Momentum 


As indicated in previous discussion, stored angular momentum is not 
an attractive source of restoring moment for the application treated in 
this report. 

As an exercise, consider the use of inertia bleeding from a wheel 
movmted parallel to the pitch axis, which is the only possible axis, given 
a geotropic pointing axis. The available angular momentum of an inertia 
wheel is, assuming maximum efficient use of materials, 

H = r (16) 

w p 

where r = radius of the wheel 
= mass of the wheel 
a = allowable material stress 
p * material density 

The average restoring moment is 


\ = H/t^ (17) 

where t^^j is the mission duration. The ratio of the restoring moments 
available from an inertia wheel and from expelled mass is, from Equations 
(9), (16) and (17) 

\ (18 

" ^sp / 
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For graphite composite, which is probably the best available material, 


a = 5x10® n/m^ 
p =1522 Kg/m® 
and \/o/p = 573.2 m/sec 

Thus, using G = 9.31 m/sec® and I„_ => 200 sec 


(19) 


In order for the inertia wheel to be competitive with expelled mass, 
it is seen that the radius of the wheel, r, should be of the same order as 
the diameter of the reflector, which is considered to be impractical if the 
required value of r is greater than the radius of the launch vehicle, as 
previously discussed. This result applies only to the average moment and 
not to dynamic variations. In the case of the latter, the momentum wheel 
is greatly superior to expelled mass because it does not run down. 

For moments whose axes are noraal to the axis of the wheel (roll and 
yaw axes) the inertia wheel acts as a g3rrostabilizer . The equations for 
motion about the roll and yaw axes are, in matrix form. 


- jA-l 

. -Hp I ( ^ ) ( ( 


( 20 ) 


where 


X = roll axis 
z = yaw axis 
(j) = roll angle 
i|< = yaw angle 
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* nioment of inertia about roll axis 
= moment of inertia about yaw axis 
p = d/dt 

H = angular momentum of inertia wheel. 


For the case of a steady rolling moment, Equation (20) may be 
solved to give 



\|; - ilssL 

H 


( 21 ) 


Both Ig and Mjj increase rapidly with vehicle size, while H is fixed 
by dimensional constraints. Thus, gyrostabilization becomes increasingly 
less attractive as the size of the vehicle is increased. Consider, for 
example, that is equal to the maximum moment produced by solar radia- 
tion pressure on a reflector with f/D =1. The moment, obtained from 
Equation (10) with Pg * .75x10”^ N/m^, is 

Mq - 1.23 X 10“^ (22) 

The angular momentum for the wheel, obtained from Eqiaation (16) with 
r = 2m, and V o/p = 573.2 m/sec, is 

H - 1146 m^ (23) 

An approximate value of the polar moment of inertia of the vehicle is 

h = (24) 
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where iHj. is the mass of the reflector. Using Equations (22), (23), and 
(24), the roll angle ^ computed from Equation (21) is: 

<(» = 1.17 X 10“^^ m (25) 

^ w 

Since 

d 2 m' (26) 

'"'Aw 

where m* is the mass per unit area. Equation (23) can be written 
w 

<f> - 1.49 X 10”^^ (27) 

K)' 

If the mass of the inertia wheel is permitted to be ten percent of 
the mass of the reflector, then 

(P = 1.49 X 10“^2 D (28) 

m' 

Both D and ra' may vary over wide ranges as indicated in Table 1. The 
most favorable values of D and m' are D = 100 m and m’ =« 1 kg/m^. In this 
case, 

» 1.49 X 10”^ radian (29) 

The least favorable values of D and m’ are D = 1000 m and m' » .01 kg/m^. 
In this case 

(|) = .149 radian (30) 
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The value given by Equation (29) is less than the smallest pointing error 
tolerance listed in Table 1, while the value given by Equation (30) is 
greater than the largest tolerance. Note that these values are instanta- 
neous values. 

The yaw angle due to steady rolling moment, obtained from Equation 
(8), using Equations (22), (23) and (24) is 

i|» ■ 1.37 X 10”^*^ ^ ^ (31) 

m* 

w 

Assuming t = 10 years = 3.15 x 10® sec, and using the favorable values 

D == 100 m and m* =0.1 kg/m^, produces 
w 

^ = 43.16 radians (32) 

which indicates the inability of the gyrostabilizer to resist steady 
moments for long time periods. 

In summary, the prescribed spectrum of mission parameters does not 
appear to favor a control system based on stored angular momentum. 
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4,3 Gravity Gradient 


Gravity gradient Is a strong candidate for control of the vehicle 
considered In this report, because the pointing axis Is geotroplc and 
because the diameter of the reflector Is large. Large diameter favors 
a gravity gradient control system because the control moment Is pro- 
portional to the fourth power of diameter, which is a higher power than 
that for any other candidate system. 

A detailed study of gravity gradient control for the present applica- 
tion is contained in Section 5. 

4.4 Aerodynamic Pressure 

Aerodynamic pressure is a candidate for control of the vehicle at 
altitudes below 1000 km. Due to geotroplc pointing, the flow velocity is 
steady relative to the vehicle, and aerodynamic control surfaces similar 
to those of a conventional airplane can be used. Pitch and yaw can be 
controlled with passive (fixed) control surfaces but roll cannot. In any 
event, the addition of active control greatly increases the effectiveness. 
For example, active trim control can reduce the steady pitch error to zero. 

If the altitude is low enough, pitching moment due to aerodynamic 
drag will be the only important component of overturning moment and for 
this case it is relatively easy to size the aerodynamic surfaces required 
to control pitch. The overturning moment for a smooth reflector is evalu- 
ated from Equation (2) of Section 3.2 and the accompanying table of pitch- 
ing coefficients vs. f/D. Selecting f/D =1.0, the pitching moment about 
the center of gravity is 
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.00277 Pd Ay • D 


(33) 


M 


where is the stagnation pressure and Ay is the area of the reflector. 
Consider a control surface of area A^, placed at a distance D/2 behind 
the center of gravity with an angle of attack a. From Equation (1) of 
Section 3.2 the control moment is 


2 

- sin^a cosa p^ A^ • D/2 (34) 

Equating to gives 

2 

sin a cosa = .00554 Ay/A^ (35) 

The maximum value of sin^a cosa is .385 at a = 54.7®. Thus, 

the minimum value of A^/A_ is 

c r 

A,/Ay - .00554/. 385 = .0144 (36) 

A somewhat larger control surface should be used to provide a margin 
of safety and also to reduce the drag of the control surface, which is 
proportional to tana. Note also, from data presented in Section 3.2, 
that the size of the overturning moment is a strong function of f/D. 

The control stirfaces required for roll and yaw are much smaller (at 
low altitudes) because the roll and yaw moments are much smaller. 
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4.5 Solar Radiation Pressure 


At altitudes above 700 lOn, the maximum overturning moment due to 
solar radiation pressure exceeds the steady pitching moment due to aero- 
dynamic drag, which suggests that solar radiation pressure may be used 
for control at such altitudes. Furthermore, in deep space, and also in 
geosynchronous orbit, solar radiation pressure is the largest of the 
environmental effects that are available for vehicle control, exceeding 
even gravity gradient. 

Solar radiation pressure is the basis of solar sail propulsion 
which has been extensively studied, References 1, 5, 6, and 7, and which 
is competitive with solar electric propulsion for high energy inter- 
planetary missions of long duration, Reference 4. 

For the present application, control moments can be provided by 
orientable reflecting panels attached to the periphery of the reflector. 
Formulas for control moments as functions of panel orientation and the 
Sun direction are published in Reference 8. Since the orientation of 
the vehicle relative to the direction of the Sun changes continuously 
during each orbit, the control surface orientations must change continu- 
ously and they must, furthermore, be actively controlled. 

At altitudes greater than 1000 Km, solar radiation pressure is the 
dominant source of overturning moment for the present application, and 
it is relatively easy to size the reflecting panels in this case. The 
overturning moment is maximum for a Sun angle of 45®, as discussed in 
Section 3.3. For f/D = 1.0, the maximum overturning moment on a smooth 
control panel is 
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- .0417 PgA^D 


(37) 


where Pg is the solar radiation pressure for normal Incidence, and 
is the area of the panel. Assume that two control siurfaces are located 
at a distance D/2 from the center of gravity and oriented as shown in 
the following sketch. 


control surface 
at zero pitch 





control surface 
at 90® pitch 



The force on the left-hand control surface is 

F =. PgAg • cos^(45*) = -i- PgA^, (38) 

and the force on the right-hand surface is zero. The resulting control 
moment is 

“c “ (39) 

Equating to gives the required size of the control surface 

A^ = .167 A^ (40) 
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A variable control moment Is produced by pitching the axes of the 
control surfaces. Since at least three control surfaces are required, the 
total control surface area equals at least one-half of the reflector area 
for f/D = 1.0. Smaller control surfaces can be achieved by placing them 
at a greater distance from the reflector, at an added cost in weight 
and complexity. Typically, the mass of a control panel, using available 
materials, is about .01 kg/m^. 

Very much smaller control surfaces are required if the reflector 
is constructed from a network of wires which produce very little net 
overturning moment due to radiation pressure. 
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4.6 Interaction with the Earth’s Magnetic Field 


If a current, I, flows through a circular loop of wire whose plane 
is parallel to a magnetic field of intensity, B, the wire experiences a 
torque 

M - I B (41) 

4 

about an axis perpendicular to the magnetic field and tangent to the 
current loop. The use of this principle for the control of orbiting 
spacecraft by interaction with the Earth’s magnetic field has been 
explored in Reference 3 and is the basis of control for a proposed low 
frequency orbiting radio telescope (LOFT), Reference 9, In that ap- 
plication the purpose of control was to produce useful scanning rates 
for a very large (1500 m) reflector. Several current loops were pro- 
vided, of which one was the rim of the reflector, and the others 
included meridional segments. Control laws were fairly complex. They 
required sensors to detect the magnetic field and computers to select 
the current loops. The weight of the control system was less than that 
of competing systems. 

Magnetic interaction has limited utility for the present applica- 
tion because it is not possible to produce a torque about an axis which 
coincides in direction with the magnetic field. In the case of an orbit 
in the magnetic equatorial plane, this axis would be the pitch axis, and 
in the case of a polar orbit it would be an axis in the longitudinal 
(roll-yaw) plane which varies with orbital latitude. Thus, to be useful, 
magnetic interaction must be supplemented by some other source of moment. 
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In order to compare the mass of a magnetic-interaction control 
system with that of other systems, the current, I, in Equation (41) should 
be expressed in terms of the mass of the control system. The mass of the 
system is concentrated in two components, a power source, probably utiliz- 
ing solar cells, and the conductors which carry the current. Assuming 
that the mass of the power source is proportional to its rated power, it 
can easily be shown that, for an optimized design, the rated current is 

D m* 

I - (42) 

8VC R’ p 

where D = diameter of the loop 

m’ = control system mass divided by reflector area 
c 

C = power source mass divided by power 
R* = resistivity of conductor material 
p = density of conductor material 

Upon inserting this value of I into Equation (39), it is seen that 

the maximum available control moment is proportional to m' D^, In this 

c 

respect a magnetic interaction control system is similar to most of the 
previously examined control systems including expelled mass, aerodynamic 
control , and control by solar radiation panels , 

In order to evaluate the mass of a magnetic interaction control 
system relative to other types, assume the following conditions 

B = 10"^ webers/m^ (typical for an altitude of 6000 Km) 

C = ,05 kg/watt 

R' = 2,83 X 10~® ohm-meter ) 

/ values for commercial 

P = 2700 Kg/m^ ( aluminum wire 
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and compute the maximum available torque. For these conditions 


M 


IT m' B ' D m^ B 
£_ ^ C 

32 C R' p 8VC R’ p 


64 X 10**5 D m^. 


( 43 ) 


For comparison, the control moment for a solar radiation control 
panel, given by Equation (39) of Section 4,5, is 

M = i. Pg A^, D (44) 

Assuming pg = .75 x 10“^ N/m^ and A^/m^ = 100 m^/kg, gives 

M = 18.75 X 10“^ D m^, (45) 

which is seen to be less than that for the magnetic interaction control 
system. Note, however, that the Earth’s magnetic field decays as 1/R^, 
so that magnetic field interaction is not competitive with solar radia- 
tion control panels at a geosynchronous altitude. 
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SECTION 5 

ANALYSIS OF GRAVITY GRADIENT CONTROL SYSTEMS 

A gravity gradient control system is particularly suitable for the 
mission considered in this report because the pointing axis is geotropic 
and because the diameter of the reflector is large. Accordingly, more 
space is assigned to gravity gradient control than to the other types 
discussed in the previous section. 

The discussion begins with passive gravity gradient control of a 
rigid body and proceeds to consideration of a control system which in- 
cludes flexible structural elements and, possibly, active control ele- 
ments. 


5,1 Stability and Pointing Error for an Uncontrolled Rigid Vehicle 


The complete small motion equations for the rotations of a rigid 
body in circular orbit are derived in Appendix D. The assumed equilib- 
rivim state is a circular orbit with the pointing axis directed toward 
the center of the Earth, and with zero angular velocity about the point- 
ing axis. This state corresponds exactly to the desired orientation of 
a geo tropic reflector. 

The conditions for stability of the desired orientation, as derived 
in Appendix D, are 


where 


^zz ^ ^xx ^ ^yy 


(46) 


zz 



, etc. 


(47) 
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and X = distance from eg in forward direction along flight path 
y = distance from eg to right, looking forward 
z = distance from eg downward toward Earth. 

For a vehicle which is symmetrical about the pointing axis, 

= lyy, 1^2 depends on the focal ratio, f/D, and on the mass placed 
at the feed. The following table shows the approximate minimum value 
of the feed mass required to satisfy the stability conditions. 

f/D °^feed 
“refl 


0.5 .50 

1.0 .07 

2.0 .015 


If the stability conditions are not satisfied, then some means of 
stability augmentation must be supplied. The simplest, but not neces- 
sarily the lightest method, is to add mass to a point on the pointing 
axis tha*’ is far from the center of gravity. Alternatively, some other 
form of stability augmentation may be used, employing any of the physical 
principles described in Section 4. 

If the stability conditions are satisfied, the angular errors may 
be found from Equation (2) of Appendix D, Since the cross-products of 
inertia are zero in the assumed load-free equilibrium state, the rota- 
tions due to small steady loads are determined by the following simple 
uncoupled formulas 


(^2z ~ ^yy) 
M 

(^zz ~ ^xx) 


(48) 


(49) 
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where n is the orbital angxilar velocity (2tt divided by the orbital 
period). Values of vs. altitude are given in the following table 
for the range of altitudes considered in this study. 


Altitude (Km) 

0 

500 

1000 

2000 

5000 

10,000 

20,000 

geosyn- 

chronous 

JJ^[(rad/sec)^xlO^] 

1.541 

1.228 

.995 

.679 

.271 

.091 

.022 

.0053 


Note from Equations (48) to (50) that it is not merely sufficient to 
satisfy the stability conditions, but also necessary to satisfy them by 
some finite amount in order to achieve small angular errors under steady 
load. Consider, for example, the steady pitching moment due to aero- 
dynamic drag, given by Equation (2) of Section 3.2, and also express 

^zz ~ ^xx 

■ =1 -i"' 

where is a dimensionless coefficient and m* is the mass per unit area 
of the reflector. For this case the steady state pitching error is 


- ^ni P d 

m' D 


(52) 


As an example, assume the following conditions 
500 Km altitude 

f/D = 1.0, and a smooth reflector 
Feed mass = 1/2 reflector mass 


Massless feed supports 


For these conditions 


= .00277 

= 20x10"^ (N/m^) 
q = .2502 

» 1.288x10"^ (rad/sec)2 
Izz/^xx = "^*98 


so that 

00 - (radians) (53) 

m D 


From the range of parameters given in Table 1, m*D ranges from 
1*0 to 1000, and the permitted tolerance on 0 ranges from •01 to 10"^# 
Thus, only for m’D > 57,3 is the steady state pointing error in pitch 
less than the maximum permitted tolerance, .01 radian. It may be con- 
cluded that, in the example considered, passive gravity gradient sta- 
bilization is adequate under some, but not all, design conditions. 

For time dependent loading conditions, the mass and Coriolis damping 
terms derived in Appendix D must also be taken into account. Dynamic 
analysis of the pitch degree of freedom is easier than roll and yaw be- 
cause its full dynamic equation is simply 

9 . My (54) 


If My is harmonically varying at frequency u, the harmonic response 
0 is given by 
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where 0^ is the static response to a load of magnitude My, and 




U, + I. 


( 56 ) 


ZZ XX 


For the example cited above, » 4.98, and o)„ - 1.413n. 

2 2 JUv y 

The amplification factor for excitation with frequency equal to the 
second harmonic of the orbital frequency is 



= -.997 


(57) 


The second harmonic of the orbital frequency is important because it 
is the dominant harmonic of solar radiation pressure. For example, if the 
orbit of the vehicle is in the ecliptic plane, then Equations (6) and (7) 
of Section 3.3 give the following equation for the second harmonic co- 
efficient of pitching moment. 

«62 ■ IS. • f Ps (58) 

The second harmonic response to a pitching moment of this magnitude 
is 

Q = S Ps . ^ 

2 3C^ m'D (59) 


The assumed value of p^ is ,75x10"^ N/m^, and, for a smooth reflector 
surface with f/D=l, Kp = ,0417. Using the parameters for the previous 
example 


02 


0.322 

m'D 


(60) 
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which has the same order of magnitude as the steady state error due to 
aerodynamic drag at 500 Km. The response to higher harmonics is much 
smaller because both the harmonic content of the pitching moment and the 
dynamic amplification factor decay rapidly with frequency. 

Nothing has been said yet about damping, and indeed no damping is 
available if the vehicle is rigid and no other means of control is pro- 
vided. Damping is necessary because otherwise transient disturbances 
persist indefinitely. Sources of transient disturbance which have been 
discussed earlier are meteoroid impact, Section 3.4, and thermal shock 
due to entry and exit from the Earth’s shadow. Section 3.1. The discus- 
sion of meteoroid Impact quoted a very small value of damping (g >■ 
7.3x10”^) required to contain the pointing error in a worst case 
analysis. 

Another matter which has not yet been discussed is the effect of 
variations in gravity gradient due to orbit ellipticity and anomalies in 
the gravitational field. The principal effect of these variations is 
to cause the stiffness coefficient represented by 3fl^ in Equation (54) 
to have a small variable part consisting of harmonics of the orbital 
frequency. The strongest harmonic of the Earth’s gravitational field 
is a second harmonic component in the meridional direction which is at 
least three orders of magnitude smaller than the steady part. Orbital 
ellipticity will produce gravitational variations which include odd 
harmonic components of the orbital frequency. 

As mentioned, an important result of such anomalies is to replace 
the stiffness coefficient in Equation (54) by 
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K = 3^2^(1 + £j^cos(Jit + <})p + e2^os(2^^t + ^ 2 ^ •••) 


(61) 


If only the first perturbation term is considered, Equation (54) becomes 
a form of the Mathieu-Hill equation which has been studied extensively 
and whose stability characteristics are well known. These studies show 
that, if ej is small, the system vd.ll be \instable if the resonant fre- 
quency given by Equation (56) lies within a small band located at all of 
the integer and half-integer harmonics of fi. For example, if =» 0,2, 
then the unstable ranges of are from .48 to .53, from .99 to 1.02, 

etc., with increasingly narrow ranges at higher harmonics. The width of 
each unstable range is approximately proportional to and since the 
expected values of the e*s are very small, all that is required to avoid 
this type of instability is to avoid nearly exact resonances with the har- 
monics and half-harmonics of the orbital frequency. This conclusion also 
holds if several non-zero values of the e^s are present simultaneously. 

A correlary of this conclusion is that it must be possible to compute the 
rigid body frequencies and some of the lower structural frequencies fairly 
accurately. At the higher harmonics of orbital frequency, very small 
amounts of internal damping are sufficient to prevent instability, even 
with exact resonances. 
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5. 2 Conceptual Design of a Gravity Gradient Control System 


It has been shown in the previous subsection that pointing errors 
can be corrected by gravity gradient stabilization, in some cases, with- 
out any explicit design modifications. This is effective for some 
vehicles, including the type considered in this report provided that it 
has a heavy feed mass. Although stability is achieved at no additional 
weight or cost, the pointing accuracy is not spectacular and damping of 
rigid body modes is not provided. 

Gravity gradient stabilization and control can be improved by the de- 
liberate addition of mechanical components to the vehicle system. Consider, 
for example, the design shown in Figure 2. The additional components are 
a small mass, called the anchor, plus a long cable and a harness consisting 
of three or more shorter cables which attach the long cable to the rim of 
the reflector. 

Each of the harness cables is in series with an elastomeric damper 
and a winch. The winches are needed only if it is desired to exercise 
active control of the pointing axis. 

Gravity gradient produces a tension in the long cable equal to 

T » 3Sl^ nig (62) 


where is the mass of the anchor (assuming that the mass of the cable 
is negligible), and z^ is the distance from the anchor to the center of 
gravity of the complete system, including reflector, feed and anchor. 
Thus, 


where m is the mass of reflector and feed 


( 63 ) 
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If the reflector is pitched through a small angle, 0, the tension 
T will produce a moment on the reflector which resists the pitch angle. 
The use of this principle to reduce pointing errors is discussed in fol- 
lowing subsections. 

A lateral force at the apex of the harness cables will tend to in- 
crease the tension in some cables and to reduce the tension in others. 
This will cause extensional motions of the elastomeric dampers (which are 
assvnned to be much less stiff than the cables) and, if the motions are 
time-dependent, damping will result. It will be shown later that both 
the rigid body modes and the elastic vibrations of the system are damped. 

Operating the winches in a differential manner will alter the junc- 
tion point of the long cable and the harness cables, thereby creating 
moments about the center of gravity of the reflector and feed. In this 
manner, the pointing axis can be deliberately altered in response to 
commands or to feedback from pointing error detectors. 

The control system shown in Figure 2 does not resist yawing moments. 
Such moments are presumably small because the reflector surface is axi- 
symmetric; also, the error tolerance for yaw is probably much larger than 
the tolerance for pointing error. The correction of yaw errors can be 
achieved by adding small masses that are connected to the rim by booms 
and dampers in the fore-aft directions. Note that booms rather than 
cables are needed because the component of gravity gradient in the fore- 
aft direction is very small (assumed zero in Appendix D) compared to 
the component in the vertical direction. 

The length of the anchor cable is likely to be of the order of 100 
times the diameter of the reflector. Based on previous work with long 
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filamentary structures, the cable should be constructed from several 
separated strands with frequent load interchanges in order to reduce 
'the probability that the cable will be severed by micrometeoroids. 
The added weight penalty to meet this threat is estimated to be low. 
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5.3 Steady State Pointing Errors 


Steady state pointing errors can easily be estimated for passive 
operation of the gravity gradient control system described in the pre- 
vious section. 

If the reflector is pitched through a small steady angle 6, the 
long cable will be translated without rotating, because the component 
of gravity gradient force in the fore-aft direction is negligibly small 
(see Equation (1) of Appendix D). As a result, the tension in the cable 
exerts a restoring pitching moment 

My = me (64) 

on the reflectors, where k is a reduction factor to account for the 
flexibility of the harness. For the present, the flexibility of the 
harness will be ignored and the reduction factor will be set equal to 
unity. The moment given by Equation (64) has an upper limit 

<Vmax “ (65) 

which occurs when 0 = D/2i. For larger angles one cable goes slack. 

An additional (positive or negative) restoring moment will be pro- 
vided by gravity gradient acting on the reflector itself, but this effect 
will be ignored because it is intended to provide a much larger restoring 
moment by adding the anchor and cable. As a result we may express the 
steady state pointing error due to aerodynamic drag by means of Equation 
(52) of the preceding section with the inertia coefficient set equal to 



where mj. is the mass of the reflector. 
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The resulting pitch angles are plotted in Figure 3 vs. the product 
Cj^m’D with altitude and focal ratio as parameters. The product m'D 
ranges from 1.0 to 1000 for the applications considered here. A prac- 
tical upper limit for is probably around 100, based on weight and 
other considerations. This is a very large Increase over the value used 
in the preceding section (.2502) for an uncontrolled vehicle with heavy 
feed mass, and ''it extends considerably the range of mission parameters 
for which passive gravity gradient stabilization is practical. 

The results for steady roll moment are similar to those for steady 
pitching moment, except that the restoring moment given by Equation (64) 
is multiplied by a factor which is approximately equal to (l+i/4Jl_). 

There is,. however, no change in the maximum moment given by Equation 
(65). 

0 

The steady pitch (or roll) error can be reduced to zero by introduc- 
ing active trim control through the winches shown in Figure 2. In effect, 
the junction point of the long cable and the harness is deflected through 
the angle, 9 relative to the reflector, rather than deflecting the point- 
ing axis of the reflector itself. As mentioned, the angle at which one of 
the harness cables goes slack is 9 = D/2Z. Thus, using Equation (52), 


fimax = U 


^m^d 

3CiA*D 


from which, using Equation (66), 



_ 2Vd 


(67) 


( 68 ) 
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This result provides the minimum value of m'Zg mg/m^ required to pro- 
duce a control moment equal to the aerodynamic pitching moment. The 
expression on the right-hand side is a function of the altitude and the 
focal ratio only. Equation (68) is plotted in Figure 4, As an example, 
if the altitude is 400 Km and f/D = 1,0, then the minimtmi required value 
of m'Zgm^/mp = 1.1. If we double this value to provide a margin of 
safety and select m* = 0.01 and Oa/nij. = 0.02, then the required cable 
length is approximately equal to 

Zg = 2 X 1,1/0,01 X .02 = 11,000 meters (69) 

independent of reflector diameter. 
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5.4 Dynamic Response 


Appendix E includes an analysis of the dynamic response in pitch of 
the gravity gradient control system described in Section 5.2. The topics 
treated in Appendix E are vibration modes, modal damping, and response to 
dynamic pitching moments. The major conclusions will be summarized here. 

The system has two pitching modes: a low frequency mode which is 

nearly a rigid body mode, and a high-frequency mode which mainly involves 
pitch of the reflector, with a small opposing motion of the anchor mass. 
Frequencies and mode shapes are plotted in Figures E-2 and E-3 vs. two 
characteristic parameters. The parameter a in these figures is directly 
related to the parameter defined in Equation (66) by 


a 



Ci 


(70) 


where I is the pitching moment of inertia of the reflector-feed system 
and m^. is the mass of the reflector. The practical range of the para- 
meter a is from about two to about two hundred. 

As mentioned in Section 5.2, damping is provided by elastomeric 
dampers in series with the harness cables. Appendix E includes an 
approximate analysis of the damping provided by the dampers to each of 
the modes. The damping for the higher mode is much larger than the damp- 
ing for the lower mode, as shown by the results tabulated in Section E.4. 
It is also shown that the damping of the low frequency mode is adequate, 
for reasonable parameter values, because the excitation of the low fre- 
quency mode is small. 
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For sinusoidal excitation, the response can be expressed as an 
equivalent static response multiplied by a dynamic amplification factor. 
The dynamic amplification factor is plotted versus frequency in Figure 
E-4 for a typical value of o and extreme values of A/Jt-. 

As discussed in Section 5.1, the most important component of sinu- 
soidal excitation is the second harmonic component of solar radiation 
pressure. Equation (59) of that section gives the second harmonic pitch- 
ing response for an uncontrolled vehicle. The same equation can be 
applied to the passive control system under discussion by using the value 
of given by Equation (66) of Section 5.3 and by replacing the dynamic 
amplification factor ( l/( l-(a)/iOy)^) in Eqxiation (59) by the value derived 
in Section E,3, From Figure E-4, a typical value of the amplification 
factor for second harmonic excitation (o)/S“2) is 1,0. Using this 
value, and previously derived parameters for solar radiation pressure, 
produces the plots of pointing error shown in Figure 5. By comparing 
Figure 5 with Figure 3, it is seen that the second harmonic response to 
solar radiation pressure is greater than the steady state response to 
aerodynamic drag for altitudes greater than 600 Km. 

The dynamic response of the system in roll and yaw has not been 
treated. The analysis is more difficult than the analysis of pitch be- 
cause the roll and yaw degrees of freedom are coupled. It is expected 
that results similar to those for pitch will be obtained, except that 
the response to second harmonic excitation may be larger due to the 
close proximity of the rigid body roll mode to the second harmonic of 
orbital frequency, indicated in Figure 1 of Appendix D, 
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SECTION 6 
ACTIVE CONTROL 


The emphasis in this report has been on the physical means uSed to 
apply control moments, rather than on control laws. Active control of a 
gravity gradient control system was discussed briefly in Section 5.2 and, 
in general , it has been mentioned repeatedly that steady state pointing 
error could be eliminated by trim control. This subject is discussed in 
some detail in Appendix F, where it is shown that integral feedback (feed- 
back which is asymtotic to a constant times l/u as the frequency, u, 
approaches zero) can entirely eliminate steady state error due to external 
(environmental) moments, and can generally reduce errors due to low fre- 
quency excitation. 

The discussion of Appendix F is largely concerned with stability. A 
practical design approach is described which yields good stability and 
good error suppression characteristics. 


58 


SECTION 7 

CONCLUDING REMARKS 

Of the sources of control moment that have been reviewed, only stored 
angular momentum and interaction with the Earth's magnetic field appear to 
be inappropriate for the mission spectrum considered in this report. All 
of the others have some potential application, depending on altitude and 
vehicle parameters. 

Due to the long mission duration, expelled mass using chemical pro- 
pulsion results in a relatively heavy control system. Electric propulsion, 
on the other hand, is a strong candidate at all altitudes. 

Aerodynamic pressure is a strong candidate at low altitudes and solar 
radiation pressure is a strong candidate at high altitudes. Gravity 
gradient is probably the strongest candidate, except at very low and at 
very high altitudes, and it is, furthermore, the only source of control 
moment that increases more rapidly with vehicle size than the environment- 
al sources of overturning moment. The requirement for geotropic pointing 
makes gravity gradient stabilization particularly appropriate for the 
mission spectrum that has been studied. 

An important aspect of control that has not been discussed to any 
extent is orbital maintenance, i.e., the application of a force at the 
center of gravity of the vehicle to maintain orbital altitude and elliptic- 
ity within acceptable limits. The only available sources of such a force 
are expelled mass and solar radiation pressure. 
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Table 1 


Mission Parameters 


Reflector Diameters: 
Circular Orbits; 

f/D; 

Physical Pointing Accuracy 
Maximum Error: 

Mass per unit area: 

Surface Solidity; 

Massc j/Mass j--, 

feed' reflector 

Mission Duration; 


100 to 1000 meters 

200 to 20,000 kilometers 
altitude and geosynchronous 

0,5 to 2.0 

lp~^ to 10~^ radians 

10”^ to 1 kg/m^ 

.01 to 1,0 
.01 to 1.0 
10 years 
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Table 2 


Summary of Control Moment Sources 


Physical Effect 

Practical 
Altitude 
range (km) 

Practical 
size range 
(m) 

Passive 

Control 

Possible 

Active 

Control 

Practical 

Expelled Mass 

Any 

Any 

Mo 

Yes 

Inertia Wheel 

Any 

< 5 

Mb 

Yes 

Gyros tab ilizer 

Any 

< 5 

Yes 

No 

Gravity Gradient 

400-20,000 

> 10 

Yes 

Yes 

Aerodynamic Pressure 

< 1000 

Any 

Yes 

Yes 

Solar Radiation Pressure 

> 600 

Any 

Mo 

Yes 

Earth's Magnetic Field 

400-20,000 

Any 

Mo 

Yes 
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Table 3 


Steady Restoring Moments for Various Control Systems 


PHSYCIAL EFFECT 



FORMULA 

DEPENDENCE ON ALTITUDE 

Expelled Mass 

M 


k G m, D 

None 



^m 

Inertia Wheel 

M 

3S 


None 

Gyros tab ilizer 

M 

a 


None 

Gravity Gradient 

M 

- 

3k 9 

RJ 

Aerodynamic Pressure 

M 


f Pd 

- o-kh 

Pd ® 

Solar Radiation Pressure 

M 

- 

‘‘s T "s 

None 

Magnetic Field 

M 

= 

-f I B 

4 

B 

r3 
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Figure 1. 
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Figure 3, Steady Pitch Angle Due to Aerodynamic Drag with 
Gravity Gradient Stabilization, 
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Appendix A 

Effect on Pointing Error Due to Thermally Induced 
_ Length Changes in the Feed Supports 


In the sketch shown below, the length of the left support is i(l + (Ae/ 2)) 
and the length of the right support is i(l - (Ae/2)) where l is the undistorted 
length and 

Ae - oAT (A-1) 

In Equation (A-1), a is the coefficient of thermal expansion and AT is 
the difference in the temperatures of the two support columns. 


— |— 



D 


A simple geometrical analysis shows that, for e « 1, the error in 
pointing accuracy due to differential thermal strain is 


9 




(A-2) 


2,/D is an easily computed fimction of the fociil ratio f/D. The following 
table lists e/6 for three values of f/D 
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f/D 


Ae/0 


0.5 1.302 

1.0 .887 

2.0 .485 

Thus, if the allowable error (due to this source) is lO”^ rad, the 
allowable differential strain for f/D =* 1,0 is ,887x10**^, Note that an 
increase of f/D decreases the allowable differential strain. 
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Appendix B 


Overturning Moment Due to Aerodynamic Drag 


The velocity vector is normal to the pointing axis as shown below 



The reflector is assumed to be a segment of a sphere and to be smooth. 
According to Newtonian theory, the aerodynamic presstire is normal to the • 
surface and proportional to cos^ where y is the angle between the 
normal to the surface and the velocity vector. The component of pressure 
in the direction of the velocity vector is equal to cosy times the normal 
pressure. Thus, the total drag force 

^d "/'Pd Yds (B-1) 

where the pressure for normal incidence, p^j »(l/^p 

The value of Cjj, the pressure coefficient for normal incidence, is ap- 
proximately equal to 2.5. The surface integral is taken over the forward 
half of the reflector. 
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The drag per unit area of the reflector is 




/ cos'^Y d S 


(B-2) 


Since pressure acts normal to the surface, the resultant aerodynamic 
force acts at the center of the sphere and may be resolved into lift and 
drag components* Only the drag component contributes to the overturning 
moment. The moment arm for the overtiirning moment is the distance from 
the center of the sphere to the center of gravity which will conservatively 
be assumed equal to the radius of the sphere, 2f. Thus, the overturning 
moment due to drag is 

Md = 2f . (B-3) 

and the moment coefficient is 

^ ^ M / cos3 V d S CB-4) 

irPjjD 

The surface integral is evaluated as follows. In terms of angles 0 and 
measured along the meridian and around the azimuth respectively, 

cosy = sln6coS({) (B-5) 


and the surface area increment 


d S » (2fr sineded^, 


(B-6) 


Thus 


/ cos^ yds 


0JJ ir/2 

4f^ I J sin^0cos'^(j)d0d({. 


o — ir/2 


— (120JJ - 8sin20jj + sin40p) 


(B-7) 


f2 

6 
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and 


(B-8) 


C 


m 



(120J3 - 8sin29p + sin40jj) 


where 


sin0jj 


D 

4f 


(B-9) 


Cjjj is evaluated below for three values of f/D. The drag coefficient 


Cd 


4F. 




2 f ^ 


(B-10) 


is also tabulated. 


f/D 

Cm 

Cd 

0.5 

.0117 

.0117 

1.0 

,00277 

.00138 

2.0 

.000672 

.000168 


The strong dependence on f/D is noted. 
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Appendix C 


Estimation of Pointing Error Due to Micrometeoroids 


A micrometeoroid which strikes the surface of the reflector will 
either penetrate it, or it will not. For particles which penetrate the 
surface, only a part of their momentum is transferred to the surface. 

It will be assumed that the maximum momentum transfer is equal to that 
of a particle which just penetrates the surface. This assumption will 
not cause Important error because the momentum transfer for larger parti- 
cles increases at a rate no greater than the diameter of the particle 
while the flux density decreases with about the third power of the par- 
ticle diameter, for particles with masses near lO"^ kg. 

The diameter of a particle which will just penetrate a sheet is 
approximately equal to one-third the thickness of the sheet (Reference 2, 
page 7), Thus, 

(C-1) 


d.. = t 


3 s 


The thickness of the sheet may be expressed as 


_ m' 


(C-2) 


where m' is the mass per unit area and Pg is the density of the material. 
The mass of the particle is, using Equations (C-1) and (C-2) 


m. 


— P (d 
6 P P'^ 


Assuming 




(C-3) 


O 

500 kg/m^ (a standard assumption) 


Pg = 2000 kg/m^ 

2 

m' =1 kg/m (the maximum range of this parameter) 
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then 


= 5x10"^ m 

dp = 1.67xl0~^ m 

fflp = 1.212x10“^ kg 

The standard velocity assumed for micrometeoroids is 30 km/ sec 
(Reference 2, page 6) so that momentum of a particle that will just 
penetrate the sheet is 

mp Vp =» 1.212 X 10“^ X 30x10^ = 3.64 x 10~^ Newton-sec 

The magnitude of the angular momentum imparted to the vehicle by 
a single micrometeoroid is 

- Xi m Vp . kp (C-4) 

where 

x^ = moment arm of impact point 
mp “ mass of particle 
Vp » velocity of particle 

kp » is a factor to account for glancing impact 

It is not possible to make a simple comparison between the angtilar 
impulse imparted by micrometeoroids and the moments due to other environ- 
mental effects. We can, however, estimate the maximum pointing error 
caused by a single impact if we make some very simple assumptions about 
the control system, namely that it acts like a simple spring to produce, 
in conjunction with the pitching inertia, I, a control frequency, oi^. 

It can then be shown that the maximum angular displacement due to a single 
impact is 
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(C-5) 


e(i) . Ji<i2 

max lo)^ 


Now I can be expressed as 

I - C. -L m' 
i 4 


(C-6) 


where is a factor that depends on shape. A typical value for is 
0.15. Substituting from Equations (C-4) and (C-6) into Equation (C-5), 
and assuming = T)l2/T 


.( 1 ) = 


®p Vp kp 


max 


Cl d3 0,^ 


(C-7) 


Note immediately that, for constant u^, 9^^^ decreases as 1/D^, so 
that the effect is largest for small vehicles.* It has already been 
shown (Equation (C-3)), that mp increases as the cube of m’ so that 
6 ^^^ increases as the square of m'.- Using the largest value of m' 

(1 kg/m'^) and the smallest value of D (100 m) , to be considered in the 
present investigation 


g(l) = 3.6 X IQ-5 

.15 X -^x 1 X (100)2 0)^. 

- 1.08 X 10"^°/o)c 


(C-8) 


*It is shown, later on, that the effect of multiple impacts increases as 
the square root of the number of impacts, which value is proportional to 
D. Thus, for multiple impacts, decreases as l/D^, 
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The number of impacts in a given time period can be estimated from the 
mean flux density, fp, of particles greater than a given size and the size 
of the vehicle. According to current estimates, the mean flux density of 
particles with mass greater than 10 kg is about 10 particles per square 
meter per second. Thus, the rate of particles of this size intersecting a 
reflector with 100 m diameter is approximately 

N* = _L d 2 f -X X (100)2 iq-7 , ^ 7 g 5 ^ jq- 3 particles/sec (C-9) 

p 4 *^4 

This amounts to one impact every 21.2 minutes which is of the same 
order as the period of pitching oscillation. 

The effects of particle impacts on pointing error are added statis- 
tically, because their directions and impact points are randomly distrib- 
uted. According to statistical theory, the combined effect of Np random 
impacts is equivalent to /Np simultaneous impacts at a point of mean radius, 
D/2/T. Thus, selecting in Equation (C-1) to be D/2/21 and kp to be 
/N^, the equivalent angular momentvnn is 

H = D mp Vp /N^ / 2/2 (C-10) 


and the maximum pointing error is 

”max 



(C-ll) 


The effect of smaller particles has not been considered in this cal- 
culation. Their effect may be estimated by assuming that the effective 
particle mass is 


m 


e 




/ “i 

o 



(C-12) 
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and that particle flux is inversely proportional to particle mass. Thus, 


where 


so that 


and 



(C-13) 


(C-14) 


(C-15) 


(C-16) 


Thus, the net effect of smaller particles is just equal to the net effect 
of larger particles. Their effect is included by doubling the value of H 
in Equation (C-10). 

It remains to select an appropriate time period to use in computing Np. 
If the control system has no damping, then the appropriate time is the use- 
ful life of the vehicle (assumed to be 10 years). In this case 


Np = • t - .785 X 10"^ X 3.15 X 10® « 2.47 x 10^ (C-17) 

2/N^ = 994.0 

and 

®max “ ^ = 1.07 x 10"^ (C-18) 

The smallest control frequency, that can be considered to be practical 
is the orbital frequency which is about 10~® rad/sec for altitudes up to 
1000 km. Thus 
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®max ” ^ = 1.07 x 10”^ radian 


which is larger than the smallest allowable pointing error (10~^ radian) 
given in Table 1, 

This result shows that the vehicle pitch mode should be damped if 
all the worst case assumptions are correct. If it is damped, the appropri- 
ate time period to use in computing Np is the time constant of the damped 
oscillations which is 

t - ^ (C-19) 

where g * 2c/c^ = 1/Q is the damping coefficient. For this case, assiun- 
ing = 10 rad/ sec, 

N = N* • t = .785 X 10”^ X 2/10"^ g = 1.57/g (C-20) 

^ P 


and 


max 


2/N^ 


3(1) 

max 


y- 


X 1.08 X 10~^ 


2.7 X 10 


-7 


g 


(C-21) 


The amount of damping required to reduce 0__„ to the minimum allow- 

max 


able value is 


(2.7 X 10"^\2 
\ 10-5 ; 


7.29 X 10"^ 


(C-22) 


which is a very small amount. 
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Appendix D 


Complete Small Motion Equations for the Rotations of a 
Rigid Body in Circular Orbit 


The equations of motion will be written with reference to a vehicle 

coordinate system whose origin is at the center of gravity of the vehicle 

and whose axes are right-handed and oriented as follows; 

X axis: pointing forward (in the orbital direction) 

y axis: pointing to the right (as sensed by an observer 

standing on the Earth, directly under the vehicle 
and looking forward) 

z axis: pointing down (toward the center of the Earth) 

The rotations about these axes are (using the right-hand rxile) 

About x-axis; <|) (roll) 

About y-axis: 6 (pitch) 

About z-axis: i(» (yaw) 


The equations of motion are derived by adding the effect of gravity 
gradient to the complete dynamic equations of a rigid body rotating about 
the negative 'y-axis with constant angular velocity, Q (the orbital fre- 
quency in rad/ sec). 

The steady inertia force density vector, including both gravity 
gradient and centrifugal force terms, is 



(D-1) 
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where R = distance from the center of gravity to the center of the Earth 
The force component, f^^, is extremely small and will be neglected. 


The complete equations of motion for pitch, roll and yaw may be 
written in the form 



(D-2) 


where p = d/dt and 


[I] 


[B] = 2n 


^yy ^zz 


-I 


xy 


-I 


xz 


”^xy 

1 ^x 

^ZZ 1 


”^xz 

T 

1 ^ ^ Vj 

L“ ! 

! 

“lyy 


j 

>> 1 
M 
1 

0 1 

~]^y_ 


Lv 1 

-^y 1 

0 



(D-3) 


(D-4) 


[K] =■ 


- lyy) 

3I*y 

"*^^xz 

3I*y 

3(Izz ■* ^xx^ 

1 ^^yz 

”^^xz 1 

^^yz 

^xx ~ ^yy 


(D-5) 



-41 


yz 


31 


xz 


xy 


+ (Terms due to other 

environmental effects) 


(D-6) 
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and where the inertia integrals 


^ab “ / » a, b =■ X, y, z (D-7) 

V 

This definition does not conform to standard definitions of moments 
and products of inertia which are, unfortvmately , not well suited for 
representing the Coriolis damping matrix. 

The inertia matrix [I] and the Coriolis damping matrix [B] are taken 
directly from Reference 4 with an appropriate transformation of axes. The 
stiffness matrix [k] has been derived from Equation (D-1) using the method 
of Reference 4. 

If the principal axes of the body are aligned with the coordinate 
axes, it is seen from Equations (D-3), (D-4) and (D-5) that the pitch degree 

of freedom is uncoupled, and that the roll and yaw degrees of freedom are 

# 

coupled by the Coriolis damping terms, ±2£2Iyyp. 

The eqtiation for uncoupled pitching oscillations is 



from which it is seen that > ^xx ® condition of stability. 

The coupled rolling and yawing frequencies are obtained by the 
solution of 

“ ^^xx'*‘^yy^ ^^zz'^’^yy^ 

- to n [4(I^jj+Iyy) (I^x'Iyy) + (^zz'*'^yy^^^xx“^yy^ ^^yy^ (D-9) 

+ 4Ji (I^^'Iyy) (^xx~^yy^ ” ^ 
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The pitching, rolling, and yawing frequencies are plotted in Figure D-1 
as functions of ^xx^^yy* ^^Sure D-1 that the 

conditions for stability (o)^ > 0) are: 

Izz > Ixx » V 
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Figure D-1. Rigid Body Modes in Geocentric Orbit 
Stability Requirement; I, , > I > I 



Appendix E 

Dynamic Analysis of, the Pitching Moments of a 
Passive Gravity Gradient Control System 

Figure E-1 shows the dynamic degrees of freedom of the vehicle in the 
pitch plane. All of the cables are assumed to be massless and inextension- 
al. The reflector-feed system is assumed to be rigid. The analysis con- 
sists of the following parts; 

!• Derive the equations of motion^ for small motions 

2. Calculate vibration modes 

3. Calculate response to sinusoidal pitching moment 

4. Estimate modal damping due to elastomeric dampers 

E.l Equations of Motion 

Referring to Figure E-1, the kinetic energy of the system is 

% = 2 ^ (E-1) 

where m = mass of reflector-feed system 
m^ = mass of anchor 

I =» inertia in pitch of reflector-feed system. 

The potential energy of the system is, for small motions, 

Ep = 2 (il 0^ + 0^) (E-2) 

where T is the tension in the long cable, given by Equations (62) and 
(63) in the main text. Note that the static restoring moment due to 
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gravity gradient acting directly on the reflector-feed system is neglected, 
as explained in Section 5.3 of the main text. 

The translational degrees of freedan, w and w. , can be eliminated 

cl 

by using the conservation of momentum 

mw + m^^w^ = 0 (E-3) 

and the kinematic relationship 

w * Wg + Z0 + (E-4) 

After eliminating w and from Equation (E-1), a straightforward 
application of Lagrange's equations yields the following equations of 
motions. 

(I + m^)p2 + Ti I WLZ^p^ 

where p = d/dt and m - mm*/(nri“m^) • 

a a 
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E.2 Calculation of Vibration Modes 


It will be noted that the determinant of Equation (E-5) is quadratic 
in Thus, the calculation of the frequencies of free vibration requires 

the solution of a quadratic equation in the variable The fre- 

quencies are expressed by the result 

(^f “ T" - 4a ] (E-6) 

N, 

where 

ft = orbital frequency 
o = mSLl^/1 
0 = 1+1 /I ^ 

The frequency ratio, u/ft , is plotted vs. o and = 0-1, in Figure 

E“2. For a > 1, which is the range of Interest, the lower of the two fre— 
quencies corresponds closely to rigid body motion with an asymtotlc frequency 
limit, oj/ft = /3", for 0 = 1.0. For a > 1, the higher of the two frequencies 
represents a pitching mode of the reflector, with little rotation of the 
long cable. 

The relative mode shapes, obtained by solving the bottom half of Equa- 
tion (E-5) for 9^/0 are given by 

“ 30(ft/u)^^- 1 (E-7) 

where u is either of the frequencies given by Equation (E-6). The relative 

mode shapes are plotted vs. a in Figinre E-3 for =0.1 and » 0.5. 

« a 
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E, 3 Response to Sinusoidal Pitching Moment 


The response to sinusoidal pitching moment may be obtained by solving 
Equation (E-5) with p = im, where u is the frequency. The resulting 
pitch of the reflector may be expressed as follows in terms of the para- 
meters a , 6 , and . 



(E-8) 


where 9^ = My/lX , is the static response. The bracketed expression is a 
dynamic amplification factor. It is plotted vs. u/Ji in Figure E-4 for a 
typical value of a and extreme values of g. It is seen in Figure E-4 that 
the low frequency mode has a narrow peak while the high frequency mode 
has a broader peak. 

Another way to calculate sinusoidal response of the reflector is to 
first compute the response of each of the modes and then to combine the 
modal responses. The following equation describes the general procedure 


uj . n % 

i k mi - 0)'^) 

where 

Fj^ = exciting force on k*^*^ degree of freedom 
mi = generalized mass of i^^ mode 
tjji = frequency of i*^ mode 
Uj = displacement of degree of freedom 

- displacement of Uj in i*"^ mode 
(j)j^ = displacement of Uj^ in i^^ mode 


(E-9) 
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For our particular case, the exciting force is a pitching moment on 
the reflector, and there are only two modes. For this case. Equation (E-1) 
can be expressed in the following particular form 

9 = ^ I ^i (E-10) 

Iuj2 i=l,2 (wi/w)^ - 1 

where the "modal participation factor", is 


i “ 1 + 0 [i/z^ + 


The following table records for o » 10 and two values of i/l. 


i/l^ 0.1 0.5 

Pj .0012 .043 

?2 .2907 .085 


Note that an Important effect of increasing i/ig^ is to greatly Increase 
the contribution of the first mode (the rigid body mode) to the response 
of the reflector. 
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E« 4 Estimate of Modal Damping and Response to Random Excitation 


The elastomeric dampers are in series with the harness cables. An 
elastomeric damper has the property that the force in the damper is out 
of phase with the motion but is nearly independent of frequency (for fre- 
quencies below 10 Hz.), Therefore, the force-displacement relationship 
for sinusoidal motions can be represented by 


F =* k(l+ig) • u (E-12) 

The per-unit critical damping of a system consisting of the damper and a 
mass is C/C^^ = g/2. 

If the damper is very stiff compared to other elements in the system 
and if g is small, then the modal damping due to the damper can be esti- 
mated from 


/modal ^ S 


(E-13) 


where Ep is the potential energy in the mode due to all stiffness elements, 
and AEp is the potential energy in the damper 




(E-14) 


If the damper is very stiff, then the force in the damper is not 
affected by the presence of the damper. Thus, in order to estimate the 
modal damping, it is only necessary to compute the force in the damper 
and the total potential energy from the undisturbed mode shape. 

For our case, the moment resisted by the dampers is the moment act- 
ing on the reflector 

M = 0)2 I 9 (E-15) 
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so that the potential energy in the dampers is 



(0^ I 02 


(E-16) 


where is the rotational stiffness of the dampers acting together to 
resist applied moment. Equation (E-2) gives the potential energy of the 
undisturbed system. Thus, using Equation (E-2), 


/C \ , g 0)^1 

r^'Vmodal 2kjTi(l + (E-17) 


where the parameters u and 0^/6 are the frequency and amplitude ratio for 
a particular mode. 

This result can be expressed as follows in terms of the parameters a> 
0, , and the ratio, 0g/9 > given by Equation (E-7) 


ls—\ 

= J.1 

fTi> 


, (o)/0)^ 1 


2l 

modal 

i'^d; 


. 9ah^ (1 + (l/^-l^ 

(0a/e)2)J 


3 /^^\« Aj 

2 l'‘d) 


(E-18) 


The following table records values of the bracketed expressed, A^, 
referred to as the "damping amplification factor", for several values of 
a and 


a 3 3 10 10 10 

Jt/ia *01 *10 *01 -10 .50 

Mode Damping Amplification Factor, A^j 

1®^ .00234 .01422 .00012 .00088 .00155 

2^° 1.0072 1.085 1.020 1.10 1.513 
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It will be noted that the damping amplification factor for the 2^*^ 
mode is high in all cases and that the damping amplification factor for 
the first (rigid body) mode is low in all cases. Since the response to 
random excitation is inversely proportion to the square root of the damp- 
ing, the above table would indicate that the response of the rigid body 
mode to random excitation would greatly exceed that of the 2®*^ mode. 
However, as Indicated in Section E, 3, the modal participation factor for 
the rigid body mode in the response of the reflector is much smaller, in 
most cases, than the modal participation factor for the 2^*^ mode. The 
ratio of the modal participation factor to the square root of the damping 
amplification factor is a good measure of the contribution of each mode 
to random oscillations of the reflectors. This ratio is recorded below 
for a = 10, and two values of 

Si 



From the data in this table, it appears that a value of l/l^ somewhere 
between 0.1 and 0.5 would minimize reflector response to random excit- 
ation. 

It should also be pointed out that the elastomeric dampers will also 
damp the vibration modes of the cables. This matter is discussed in 
Reference 1. 
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Figure E-1. Motions of System in Pitch Plane. 
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Appendix F 

Elimination of Steady-State Error by Integral Feedback 

In aircraft, steady-state pointing error is eliminated by trim con- 
trol, which can be regarded as a very low frequency control system, often 
with the pilot in the control loop. The subject of low frequency control 
is investigated in this Appendix from the viewpoint of control system 
theory. 

Figure F— 1 shows a conventional block diagram for a control system 
interacting with a rigid body, whose rotational inertia is I. In Figure 
F-1, is an input command, and M is an external (environmental) moment 

applied to the rigid body. For no input command, it is easily shown that 
the pointing error is 

e - ^ (F-1) 

I + F(p) 

wherein p may either be regarded as frequency (p = im) or as a deri- 
vative operator. In order to achieve zero error as p -»■ 0, it is clear 
that the feedback F(p) must be asjnntotic to A/fP’ where A and n are con- 
stants and n > 0, In discrete systems, n is an integer and F(p) is the 
ratio of two polynomials, F(p) = G(p)/H(p), The control system design 
problem is to select F(p) so that e is smaller than a desired upper 
limit which may be a function of frequency, and such that the system is 
stable. 

Stability is determined by examining the roots of 

I p2 . H(p) + G(p) = 0 (F-2) 
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Elementary analysis shows that the simple forms 


F(p) » A/p 

and (F-3) 

F(p) = A(p+Pq)/p 


produce unstable roots.* The simplest form that will produce stable roots 
is 


F(p) = A(p2 + (o^gp + 0)2) 

P 


where A, and g are constants. The asymtotic form of the error for 
p 0 is 


Mp 

e Y 

Au)^o 

2 

which indicates the desirability of a large steady state gain, Awq. 
The maximum error occurs near p = im^ and has a magnitude e = M/A<UQg. 
The stability equation for this case is 


(F-5) 


3 2 2 
Ip + A p-^ + Am gp + Am_ 


(F-6) 


and the conditions for stability (roots with negative real parts) are that 
A and m^g be positive and that 


A > 



(F-7) 


2 

*For the reason that the coefficients of p and must be positive 
(not merely zero) for a cubic equation to have all stable roots. 
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Equation (F-7) indicates a lower limit for the gain, A. An upper 


limit is provided by the requirement that the system be stable in the 
presence of structural vibration modes. In order to include such modes, 
the impedance of the structure Z = 1/lp^ is replaced by 

Ip^ i Ii(p^ + SiWiP + (F-8) 

where Ij^, Wji, and are respectively the generalized mass, frequency, 
and damping of the 1*"^ mode. The loop gain of the system is 

L » F(p)/Z (F-9) 


Assuming small structural damping, the loop gain will have peaks at 
the modal frequencies, , If » tOg then, from Equation (F-4), 

F(io)j^) ~ iAoj^ near a structural peak and the i^^ peak in the loop 
gain has the approximate magnitude 




_A 

IigiU>i 


(F-10) 


In order to provide gain stability relative to every structural mode, 
it is required that L(aij^) < 1, or 

A < liSiWi (F-11) 

Combining this stability condition with the low frequency stability 
condition. Equation (F-7), gives 

< A < (F-12) 
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which can only be satisfied if 


^i“i 


g Si > 1 


(F-13) 


It is difficult in practice to ensure that this condition will be 
satisfied. The control system can be improved by adding high frequency 
roll— off which will ensure that none of the structural modes are de- 
stabilized. The roll-off must be at least second order in order to 
satisfy the physical requirement that the feedback F(p) be zero at 
infinite frequency. If such roll-off is not provided deliberately, it 
will be provided by structural flexibility and other parasitic effects, 
with potentially unpleasant consequences. With the addition of second 
order roll-off, the feedback function is 

.,2 2 . 

A(p + ajjjgp + oJq) 

“ P(l+p/a)l)(l+p/m2) 


which is the simplest possible form for practical integral feedback* 
2 

If 0 )^ « a)j 0 ) 2 » then it can be shown that the conditions for 
stability are 



< A < I (o) 1+0)2) 


(F-15) 


As a matter of practical design. A, o)q, o)i» o) 2 > and g should all 
be chosen as large as possible in order to reduce the low frequency errors. 
The values of these parameters are, however, inevitably limited by the 
presence of structural modes. 
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CHAPTER 2 


EFFECTS OF PHASE ERRORS ON ANTENNA PERFORMANCE 


by 

3ohn M. Hedgepeth and Karl Knapp 
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SECTION 1 
INTRODUCTION 

The successful performance of any structure depends largely 
on the identification of the critical or primary loads and design 
criteria on which the design is based. The future promises large 
structures which must be deployed, erected, assembled, or fabri- 
cated in space. For such structures, which will not be required 
to face the launch environment, the primary design requirements 
will be derived from the space-flight environment and will deal 
with phenomena as primary criteria which have been considered as 
only secondary in the past. The design of such genuine "space" 
structures will require a solid foundation of critical criteria. 

This chapter examines the influences of deviations of the 
surface of a large radio-frequency antenna on the antenna pattern. 
Clearly, such information is needed in order to establish require- 
ments on the accuracy of the structure of the antenna so that it 
can properly perform its mission function. 

From a mission-usefulness point of view, three characteristics 
of an antenna pattern are important. The first is the intensity 
of the radiated pattern at its center. This "on-axis gain" or 
"main-lobe gain" is important because it establishes how much power 
must be supplied to the antenna in order to achieve an acceptable 
power density at the receiving station. Alternately, if the 
antenna is used as a receiver, the main-lobe gain determines the 
strength of the received signal relative to the strength of the 
various noise sources in the antenna and other parts of the re- 
ceiver system. 

For some missions the amount of energy of the main antenna 
lobe is important as the main-lobe gain. This so-called "side- 
lobe gain" is of concern in those type of missions in which the 
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radiation of energy outside the main lobe is undesirable. Geosyn- 
chronous communication systems, for example, must have low side 
lobes in order to avoid crosstalk among the various satellite- 
Earth station links that are using the same frequency channels. 

In sensor applications, low side lobes are required in order to 
eliminate the ambiguity that arises because of off-axis "hot 
spots . " 

The third characteristic is that of "main-beam efficiency." 
Antennas used to transmit power, for example, must concentrate 
the power in the main beam so that the receiver can capture it. 

A sensing antenna would need to minimize the power gathered out- 
side of the main beam if the contrast were low. Microwave radio- 
metry is an example of such an application. 

The effect of antenna errors on the gain pattern arises pri- 
marily from nonuniformity of the phase of the radiated signal as 
it crosses the aperture plane of the antenna. This phase error 
in a reflector type of antenna is composed of errors arising from 
the antenna feed, those due to blockage, those caused at the re- 
flector surface, and those due to improper geometry of the re- 
flector and feed position. (If the reflector were a perfect 
paraboloid and the feed were located exactly at its focal point, 
there would be no errors from antenna geometry.) For directly 
radiating antennas, a similar separation of phase-error sources 
is possible: those due to inaccuracies in the signal fed to various 

antenna elements and those due to imperfect antenna geometry. In 
this paper, attention is focused on the phase error itself without 
consideration of the source; of course, the overall intent is to 
investigate these types of phase errors which might arise from 
geometrical imperfections. 

The loss of main-lobe gain in the far field due to phase im- 
perfections is generally agreed to be a straightforward function 
of an average phase nonuniformity. To be precise, the main-lobe 
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2 2 

gain is proportional to exp (-6 ) where 6 is the weighted mean 
square deviation of phase (in radians) , the mean being taken 
across the aperture with a weighting function of the illumination 
intensity. This formula is given by Ruze in his classic survey 
paper (ref. 1) . It agrees for small values of 6 with previous 
work by others, as well as the earlier work of Ruze himself. In- 
cidentally, Ruze's earlier work contains a different form for 

2 

larger values of 6 and purports to show that the main-lobe gain 
depends on the distribution of the error across the aperture in 
addition to the weighted mean square (ref. 2) . The more recent 
result, in which the areal distribution has no independent effect 
on the main-lobe gain, is used in the results herein. 

The side-lobe gain is generally agreed to be influenced by 
the distribution of the errors as well as their weighted mean. 
Ruze shows, for example, that in the case of random errors, the 
side-lobe gain is a strong function of the correlation interval 
(ref. 1) . Different types of structural phenomena will create 
different types of phase-error distribution. For example, the 
distortions due to thermal effects or maneuver loads are likely 
to have characteristic lengths that are a large fraction of the 
aperture size. Such large-scale errors are termed herein as 
"overall" errors. On the other hand, errors due to random manu- 
facturing inaccuracies, which are treated in ref. 3, generally 
have much smaller characteristic lengths. Such small-scale 
errors are termed herein as "local" errors. 

The loss in main-beam efficiency due to phase errors is 
dependent on the effects on main-lobe shape as well as main-lobe 
gain. 

This paper contains newly derived results as well as results 
taken from and developed from the literature for the effect on 
main-lobe gain and shape and side-lobe gain of various types of 
imperfection patterns. 
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SECTION 2 


BASIC ANTENNA FIELD EQUATIONS 
AND PHASE ERROR-FREE SOLUTIONS 

An estimate of the scalar diffraction pattern for a plane 
circular aperture can be obtained by using the Fraunhofer approxi- 
mation : 


EO.M - 8 cos 

*X) ‘'o 

where f(r,4>') is the magnitude of the illumination at the aperture 
and 6(r,4>’) is its phase. The coordinates in respect to a parabolic 
reflector antenna are shown in Figure 1. The above expression is 
referred to as the far— field equation, as it is a good approxima- 
tion of the antenna field only when the distance from the antenna 
is large compared with the aperture radius a and beyond the inter - 
fsrence pattern of the Fresnal zone. Reference 4 contains a com- 
plete discussion of the derivation of the Fraunhofer approximation. 

For a circular aperture, it is convenient to change the vari- 
ables so that 

2iTa . „ 

u = — r — sin 0 


and 


P 


a 


so that the above equation becomes 
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E(u,(j)) 




pdpd(|)' 


( 2 ) 


= a 


2 yi" ^ 

y rf(p,(}>')e^t<S(P,<})')+up cos 

^r\ n 


0 *^0 


If we restrict ourselves to the circularly symmetric feed il- 
luminations characteristic of this type of antenna, then f will not 
depend on (p ' so that the field equation becomes 


E(u,« - ayyf(p)e^t*<P-*'> + "'’ pdod*- (3) 

-V 

In the special case where there are no phase errors, 6=0, the 
above expression can be integrated with respect to (p' which results 
in 


E(u) 



f(p) Jq(uP) pdp 


(4) 


For a uniform illumination, the function f = 1 and the far- 
field pattern become 


E(u) 


2TTa^ 


(u) 

u 


(5) 


The equation for the more general case containing phase errors and 
tapered illuminations is more difficult to integrate. 

The gain of unblocked antenna is given by 


G(0,(j)) 


4 ^ yF(e,r)dA|^ 
|f( 9,(J)') l^dA 


( 6 ) 
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where F{0,(t)') = f (p , c}) ' ) ^ ^ ■*■ (4>-4>')] and A is the 

area of the aperture. 

In the special case, mentioned above, of a unifoinn illumina- 
tion with no phase errors, that gain is given by 







(7) 


and the maximum gain occurs on the z axis where u = 0. 





( 8 ) 


This represents a theoretical maximum for a paraboloid antenna. 

The effects of nonuniform illumination, phase errors, spillage, and 
blocking by the feed and its supports all reduce the gain achieved 
in practice. 

Feed systems produce nonuniform illumination of the reflector 
resulting in nonuniform aperture fields. Generally, these distri- 
butions peak at the center of the reflector and taper to a low 
value at the edge. For present purposes, it is assiamed that none 
of the feed illumination spills over the reflector edge. Spencer 
(ref. 5) has treated the case of distributions of the form 
fa(l-r )^, the treatment resulting in the general solution for no 
phase errors 


E(u) = 


^2 2 
TT a 


2^ pi J 


p+1 


(u) 


.P+1 




(u) 


(9) 
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where A is the Lambda function. The characteristics of these 
diffraction patterns are shown in Table I. The solution for p = 0 
is the uniform illumination case. These tabulated results indicate 
a reduction in the gain and a widening of the beam as the distri- 
bution becomes less uniform over the radius of the circular aper- 
ture. However, they also indicate that intensity of the first side 
lobe in relation to the on-axis gain is significantly reduced as p 
is increased. The gain factor, G , is the on-axis gain divided by 
the maximum on-axis gain of Eq. (8) for p = 0. 
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SECTION 3 

THE EFFECTS OF OVERALL PHASE ERRORS 

One important class of phase errors results from overall dis- 
tortions in the mechanical shape of the antenna reflector. In 
order to examine the effects of these errors, circularly symmetric 
illuminations of the form 

f(p) = (l+p)(l-p2)P (10) 

will be considered where p is an integer. Note that f(p) has been 
normalized so that 



( 11 ) 


Phase errors of the form 
5(P,4>') = g(p) cos n(j)' 

will be considered. 

The mean square error is defined to be 


T. L 


f(p)r(p,(|) ')dA 
A(p)dA 

Jk 


( 12 ) 


(13) 


One method of solution is to expand the exponential term of 
Eq. (2) for very small errors so that 
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E(u,4.) : f r'° 

•^0 /d 

x|^r + ig(p) cos n(})' - ■ ^ -■ | --- cos^ n(|) 'J dpdcl) ' (14) 

2 

where g(p) <<1 and terms higher than g (p) have been dropped. 

Let 

(0 = 4)' - <j) 

then 



By expanding the trigonometric terms in ((f)+w) and noting that the 
sine terms integrate to zero, since e^^ cos w even function, 
then 
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E(u,(t)) 


- 1 2TT-(f) . 

a^ypf(p) y 

0 -(p 

X 1^1 + g(p) cos n(f) cos no) - ^ cos 2n<p cos 2ncoJ dpdw 

(16) 


This expression can be integrated with respect to the angle o) so 
that 


E(u,(j)) = 2TTa‘ 


y'pt(p)[(i-i^) 


Jq(up) + 1 


(n+1) 


g(P) Jj^(up) 


cos n({) 


3.1 UNIFORM ILLUMINATION AND LINEAR PHASE ERRORS 

As a first case, consider uniform illumination and linear phase 
errors : 


Uniform illumination, f(p) = 1 

Linear phase errors, 6(p,4>') = 2v 6 p cos n<j)' 
Substituting these expressions into Eq. (18) results in 
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E(u,(l)) = 2A 


2/^3 

,(up)dp-6 r p jQ(up)dp 

'0 *4 


f 


+ iJ 6^ ' cos nd> f 


cos n(J> /*p J (up)dp 

J “ 

•'o 


+ (-1)^"*"^ 6^ cos 2n(}) C p^J^(up)dp 


(18) 


This expression can be evaluated by series expansions where 


? = f (up)dp = 

^ Jr^ ^ __0 (2p-hi+3)p!(p+n).' ' ' 


2p+n 


1 ” ^ 

’’a ’ f - 2^ 

*'o p=0 


(- 1 )' 


(2p4n+4)pl (p+n) 


;•( 


\2p+n 

t) 


(19) 


A program suitable for use in a TI-58 or -59 calculator is 
shown in Appendix A. This program was used to evaluate these terms. 
Note that 


1 2 

■^ = I P Jj^(up)dp = 

*{) 


u 


r 3 

\ = f P 

*{) 


U 


( 20 ) 
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Also note that 


11m j. , _ 1 

u-K) ^0^“^ 3 


lira ^ _ 1 

= 4 


and that, in general, for n > 0 the limit of each function is 
as u goes to zero. 

To examine one case in particular, let n = 2 so that 


6(P 


“ iJ 6 ^ ' 


6 p cos 24 »’ 


Then 


E(u,(j>) = 2A 


1 ^- 


6^ [Gq(u)+G^(u) cos 2 / (u) cos 2(j)] 


and 


E^(u,(}>)| = 4A^ 




6^ (Gq4G^ cos 44>)J +4 6^F2 cos^ 2<j)| 


which, after dropping terms higher than 5^, reduces to 


lE^I = 


4J, -r 8J 

-T-« 

u 




I+G^ cos 44> - 


2uF; 


2 2 


cos 2<}) \ 

/ 


(21) 


zero 


( 22 ) 


(23) 


(24) 
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Since 


/ 


F dA = A 


(25) 


the gain is given by 


G(u,(|)) = 


4ttA 

,2 


~1 / 2uF \ 

6 “;7“(Gq+G^ cos 4<j) — cos 2(j)j 


(26) 


which has been plotted in Figure 2 for/ 6"^ = 0, 0.125, and 0.250 
when 4) = 0. This expression is identical to Eg. (7) for 6 = 0. 

3.2 QUADRATIC ILLUMINATION AND LINEAR PHASE ERRORS 

As a second case, consider quadratic illumination and linear 
errors : 

2 

Quadratic illumination, f(p) = 2(l-p ) 


Linear phase errors, 6(p,4>') = V^v 6"^ p cos u<p^ 


Substituting these relations into Eq. (18) results in 


E(u,4>) S 27Ta^ r (2p-2p^) 
*{) 




2 2 

1 - I (up) + 1°'*'^ V^/ 6^ p J^(up) cos n(}) 


' , , 1 ^n+1 36^p^ T / V 

+ (-1) — 2 — “^2n 


dp 


(27) 
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If we define 


H 

n 


I 

n 



J^(up)dp 


J^(up)dp 


(28) 


Eq. (28) can be evaluated using the program in Appendix A. Then 


E(u,4)) = 4A 


|v-°o 


+ 6^ (F -H ) COS n(}) 

n n ^ 


- 6" 




<'^2n'^2n> "* 


•1 


(29) 


Since 


/ 


? , 1 2tt 

iF^ldA = /* C 4(l-2p^V)pdpd(()’ = f A 


) ^ ^ 

ff 

0 


(30) 


the gain for n = 2, after neglecting terms higher than 6"^, is 


r/ ~ 3 ttA 

G(u,<{>) = 


I (31 


- 2(F2~H2)^ cos^ 24) 


These results are shown in Figure 3 for/ 6'^ = 0, 0.125, and 0.250 
when 4> = 0 . 
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3.3 QUADRATIC ILLUMINATION AND QUADRATIC ERRORS 

A third case of interest is to examine the effect of quadratic 
errors combined with a quadratic illumination. 

Quadratic illumination, f(p) = 

Quadratic phase errors, 6(p,<{)') 

Substituting these expressions into Eq. (18) results in 


2(lV) 


= 2 /3/ 6^ p^ cos (})’ 


E(u,<()) = 2TTa 


r 2(p-p^) (1-3 6^p^) Jq(up) - 3 6^p^J^(up) cos 44> 


2 2 . 


- i 2 v3v 6 p J 2 (up) cos 2(f) | dp 


■] 


(32) 


If we define 


K = r Q 3 
n I n 


(up)dp 


(33) 


Eq. (33) can be evaluated with the program in Appendix A. Then 
for (j) = 0 


E(u,0) = 4A 


[v- =0 - 


i 2 - 6^ 3(1 




(34) 


and 


0(u,0) - ^ |l 6 I 




(35) 


These results are shown in Figure 4 for v 6 =0, 0.125, and 

0.250. 
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3.4 QUADRATIC ILLUMINATION OF A SPHERICAL REFLECTOR 


The method of solution used in the above three cases is lim- 
ited to very small values of/ 6^ in order to ensure that g(p) << 1. 
A second technique can be used which avoids this limitation by per- 
forming a numerical integration of the far-field equation. 

Consider phase errors of the form 


6(P) 



(36) 


which result from the use of a spherical reflector instead of a 
paraboloidal reflector, and where 


^ • 6'^ JJiJ 



(37) 


where F and D are the focal length and diameter of the reflector, 
respectively. This relationship applies only to the case where 
the feed is located at the focal point. The errors can be reduced 
for a given wavelength and geometry by displacing the feed along 
the z-axis. 

Substituting Eq. (36) and a quadratic illumination into Eq. 

(3) results in 

O «iy 1 / 0 / 

E(u,*) = /(P-P')e^ ‘ (p'’-l/6)^iup cos 

* 1 ) /) 

This expression can be integrated with respect to (j) ' so that 
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2 r (r. ^ ^ ^ 

E(u,(f») = Alta / (p-p )e '^jQ(up)dp 

*^0 


(39) 


which, after expressing the exponential in terms of trigonometric 
functions , results in 


E(u) = TTa^ A ^ (pV)Jq(uP) 


COS 


^ 6 V^ 5 / tJ ^ j cos 1 ^/ 5 / 7 ‘v^ 6 ^ 'j 


+ sin 


^6 45 n J 6 ^p^) sin (v^ 5 / 7 V^ 6 ^ j 
^6 / 5 / 7 >/^ sin 
- i cos ^6 y/bjiJ 6^ ^jcos 


+ i sin 



(40) 


Define two functions 


M(u) - j (p-p^)J^(up) cos /V 7 ‘)/^^p^jdp 


'I 


(< 


N(u) = A 


/ (P-P^)Jq(up) sin (e V^/^A)dp 

*^r\ 


(Al) 


so that the field expression becomes 


E(u) = A ( M cos >/ 577 ‘ / 6 ^ -N sin / 6 ^ 


+ i N cos \/ 5 T 7 / 6 - i M sin a/s/?' / 6 



(A2) 
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The gain, according to Eq. (6), is then 


G(u) = ^ (mW) (43) 

The terms M and N can be evaluated for various values of u 
by numerical integration. Appendix B includes a program suitable 
for use in a TI-58 or -59 calculator which was us ed t o obtain the 
values plotted in Figure 5 for several values ofVo^. 
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SECTION 4 


THE EFFECTS OF SMALL WAVELENGTH ERRORS 
ON THE SECONDARY PATTERNS OF ANTENNAS WITH CIRCULAR APERTURES 


To establish the relative influence of small-scale phase 
errors, consider the case of uniform illumination and linear phase 
errors previously considered in Section 2.1. Now, a solution will 
be obta ined for large n. Substituting f(p) = 1 and 6(p,(|)') = 

2 V 6^ p cos u<t> ' into Eq. (2) results in 


E(u,<})) 


27T 1 . 

-'a 


6 pcos n4»'^iup cos 


'0 ''0 


(44) 


which, after substituting a = becomes 


E(u.4.) . ^ f 

0 


r i- 2v 6^ p cos a iup cos (<))-a/n) 


pdpda 


(45) 


The intergral over a to a limit at 2mr can be replaced by the 
summation of n integrals to a limit at 2 tt. 


2 1 


E(u(|)) ~ pd p^ r 

p=l*0 


2 tt . -/^2 
* - 2v 0 p cos a 


^ g(iup cos a/n + 2Cp-l)iT/n - (|))^^ 
for large n and a< 2ir, a/n << 1. 


(46) 
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Therefore, the asymptotic approximation is 



(47) 


Since 


n 


^ V ■' iup cos 

n / j 
p=l 



1 r iup cos 
0 


(X-(|)) 


dX 


Jq(uP) 


(48) 


then 


E(u,<J>) 



" 0^2 


P) jQ(up)pdp 


(49) 


which, when evaluated, becomes 


E(u,(|)) = iTTa^ — ^ ■ - 

The resulting gain pattern is 

G(u) : ^ ^|u J|j(2/?)Jj^(u)-2/?j(|(u)Jj^(2/«^ 

(u^- 2 ?) 

The on-axis gain is 

jJ(2/?) 

G(0) S -i 


u J (2 
0 


6‘-)Jj^(u) -2 


6^ Jq (u)J^(2 


6 ^) 


(50) 


(51) 


(52) 
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and the ratio of gain to on-axis gain is 


G(u) 

G(0) 


4J 

(u^-46^)^ 


Jq(2v' 6-) f=^ 

2 =rr - 2/ 6^ Jq(u) 


J, (2/^) 


The gain pattern has been numercially evaluated for several 
of V 6^ and the results are shown in Figure 6. 


(53) 

values 
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SECTION 5 


THE EFFECTS OF PHASE ERRORS 
ON THE SECONDARY PATTERN OF A SQUARE ANTENNA ARRAY 


The far-field equation for a square antenna array with side 
of length L can be written as 


L/2 L/2 

-L/2 -L/2 


(54) 


Consider the aperture to be divided into an N x N array as ele- 
mental squares, each one of which is denoted by integers m and n 
(1,2,3...N). In the elemental square (m,n) , let 

C = F + (2in-l)C-^ 
m I 

n = n+(2n-i)c-y 
n 2 


where 


C 


JL 

2N 


and let f and 5 be constant. Then 
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N N 




^ ^^^inn^'-(27Ti/X) (L/2+C) sin 9 (cos ({) + sin (})) 

mn 


m=l n=l 


X e 


(47TiC/X)sin 6(m cos (f) + n sin (|)) 


C C 

.ffe 

-c -c 


(27ri/X) sin 9(5 cos <}*+n sin <))) 


d5dri 


(55) 


Integrating and squaring the absolute value yields the power to be 


E(6,<{)) 


• /27rC . Q , 

sinl— ^ sin 6 cos q> 


jsin ^ — sin 9 sin 


^ sin^ 9 cos <J) sin (j) 


N N N N 

EEEL 

m=l n=l p=l q=l 


f f e 
mn pq 


i(6 -6 ) 

mn pq 


X e 


4TTiC/A sin9[(m-p) cos (|)+(n-q) sin 


(56) 


5.1 UNIFORM ILLUMINATION AND NO PHASE ERRORS 
If f = f and 6^^ = 0, then 


mn 




mn 


sin|^^ sin 9 cos <J>jsin|^^^ sin 9 sin (^j| 


tt2 2 

— T sin 9 cos (}) sin 4> 


. / 21 TCN . n A . /2itCN . r. a 

sinl — ^ — sin 0 cos (plsinl — ^ — sin 9 sin <pl 

. /2ttc ; I ,\ . /2irc ; I rr 

sin I ■ sin o cos 91 sinl — sin 9 sin (pj 


(57) 
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(4cV)^ 


(58) 


|Eq(0.0) 


= L" 


Let the relative power be 


P(9,(j)) = 


E(e,d)) 

|Eq(0,0) 


(59) 


Then 


p(e,(j)) = ^ 


sin 6 cos sin 9 sin 


7T^ 2 

— r sin 6 cos (j) sin (f) 


N N N N 

ESEi: 

m=l n=l p=l q=l 


£ f ^^^mn *^pq^ (AiriC/X) sin 9 
tie e 

mn pq 


X (m-p) cos (p + (n-q) sin cp 


(60) 


For uniform illumination and 6 = 0, 

mn ' 


Po(e,(|)) = 


sin^^ sin 9 cos <j)jsin|^ sin 9 sin (j| 


2 2 

L . 2 Q A ^ A 

— ^ sin 9 cos <p sin cp 


(61) 


For 4> = 0 , tt/4 
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PqO.ttM) 


. 2/ TL . Q,| 

sin sin 0)1 

2 2 

Tl . 2 . 

:r- sin 0 

L 2A^ 


(62) 


These gain patterns are shown in Figure 7 . 

5.2 UNIFORM ILLUMINATION AND RANDOM ERRORS 

Let be a random variable and let the illumination be uni- 
form. Then, the expectation of the complex exponential is (see 
ref. 1) 


/ i(6 -6 

/ mn pq 



1 when m = p and n = q 



when m 5^ p or n q 


(63) 


where is the standard deviation of 6, assumed to be uniform. 
Note that statistical independence is also assumed. 

The corresponding expectation of the gain is 


<^P(0, ().))> = Pq(0,()))< 


-a. 


.Ld) 


N 


. /zttC . q a . /2ttC . a . a 

sxnl— ^ sin 9 cos (pi sinl-^^ ® sin (pi 


( 






sin 6 cos ()) sin ({> 


(6A) 
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The main-lobe gain is 


-a 


-o. 


GflCS) 


= P(0,0) - e 


6 ^ 1-e 


(65) 


N 


Define the additional side-lobe gain due to errors as 


S.L. 


P(9,(|)) -PQ(0,(j))GQ(6) 


( 66 ) 


For the present case, this gives 


-o. 


1-e 


N 


S.L. 


^2 

—Go - 0 

0 . 1-e 

e + r — 


sin 9 cos (f>j sin 6 sin (})j 


(■ 


sin'^ 0 cos (j) sin (J> 


Po(6.4>) 


N‘ 


(67) 

Of course this is zero for 0=0. If we look at the region 
of the first side lobe and consider C<<L, then we have 



There is a grating side lobe (for example, at ttC/X sin 0 = 3 tt/ 2 
for (f) = 0) but its magnitude is smaller than that above. 
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5.3 UNIFORM ILLUMINATION: A CHECKERBOARD 

PATTERN OF PHASE ERRORS 


6 = 6 , m + n even 

mn 


= -6, m+n odd 


We want to evaluate the suminations in Eq. (59) for uniform illumin- 
ation. First, consider a single summation 


i 


C/X)sin 0m cos (J) _ ^i6(-l) 


m=2 ,4,6 


(4tt 1C/X) sin 0m cos <j> 




I IN--L 

E' 

iif=1,3,5 


(AiTiC/X) sin 9m cos <j) 


Note that we have assumed that N is even. 

The single summation can be closed to yield 


^i6(-l)^ ^ ^ ^-2i6 ^_^^n^-(4TTiC/X)sin 0 cos ^ (N+2) (2'iriC/X) sin 0 cos <}> 


. /2N1TC . 


sin 0 cos <})) 


sin|^^^ sin 0 cos (fij 
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The double summation over m and p, after some rearrangement 


N N 


m=l p=l 


pq^ (A'TTiC/X) sin 0(m-p) cos (j) 


. /ttL . „ ^ 

sin 1^ sin D cos (pi 

= 2 /TtT; V [cos a 6 + cos(b 6 + 4 itC/A sin 9 cos d))] 

sin sin 6 cos (}.) nq nq 


where 


n even 


n odd 
n odd 


q even 
q odd 
q even 
q odd 


Now, the sum over n is 




-1) (4TTiC/X) sin 0n sin <J) 


-t sin ( ^ sin 0 sin <j) I 

- + ♦J — p,K , a ■ 

L sin sin 0 sin <p 1 


So the quadruple summation is 
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Ag = Pq( 0,<{)) sin 0 cos (J)|tan^p^ sin 0 sin (f)jtan^ 6 (76) 

This is zero for ((> = 0 and = tt/2. For other values of ip 
and for the first side lobe with L>>C 
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(77) 


sin 


[ttL 

i X 


sin 0 cos (j) 


\ . 2/ttL 

7!!M~ 


sin 6 sin 4>j 


F.S.L. 


tan 


N 


The location of the first side lobe in the octant (|) < tt/4 is 
approximately at 


sin 0 cos (j) 


3tt 

2 


So 


A 


F.S.L. 



(78) 


This is largest at <}) = tt/4 . Thus 


^F.S.L. 4 ^ 

N 

This is fairly small for reasonable 5 and N ^ 10. 

There are grating side lobes located, for example, at 


. q , IT 

sxn 0 cos (}) = y 

At this location 


(80) 


TTL 

A 


sin 0 cos (j) 


Lit _ NTT 
4C " 2 


(81) 


Since N is even, Pq( 0,<()) = 0. Thus we have a indeterminate value 
for the grating side-lobe gain. 
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The limit is 


, • 2irC , a . 
Ixm — r — sxn 0 cos <p 


/til 

^(t 


sinl-Y" sin 0 cos (|> 


2/ttC . . . 

cos sin 6 cos <p 


-)J 


= N(-l) 


N/2 


(82) 


(83) 


Therefore 


G.S.L. 


2 2 2 
= 16 sin (Ntt/ 2 tan (})) tan (tt/ 2 tan (}>) tan 6 

N^TT^ tan^ 6 


(84) 


Again, this is maximum for tan 4> = 1 and is 


G.S.L. 


^ 2 . 
— 7 * tan 0 

TT^ 


(85) 


which is not small. 


<() = 


IT 


G(9,tt/ 4) = 1^1 -sin^ 6 - tan^ 

where 
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Z = sin 0 = — (86) 

/2A /2 


The secondary pattern has been numerically evaluated and the 
results are displayed in Figure 8. 
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SECTION 6 

CONCLUDING DISCUSSION 


The effect of phase errors on main-lobe gain is summarized in 
Figure 9. The narrow band labeled "random error” applies to both 
the analysis of Ruze (ref. 1) and the square aperture results of 
this paper (Eq. 65) . The "scallop error" and "periodic checker- 
board error" curves are obtained from the large— n and square— aperture 
analyses herein (Eqs. 52 and 75). The "1-6^" approximation was 
obtained early by Spencer (ref. 5) . All the other curves differ 


little from e 



Examination of the preceding beam patterns indicates that the 
errors affect primarily the magnitude of the main lobe and only 
secondarily its shape. Therefore, the effect of errors on the 
main-lobe gain is also interpreted as the effect on main-beam ef- 
ficiency. 


Also shown in Figure 9 are the values of rms phase error that 
are associated with rms surface errors of a reflector antenna of 
A./50 and X/16. Note that X/16 is usually considered to be adequate 
accuracy for point-to-point conununication insofar as the aperture 
size can be increased to achieve the desired gain. However, if 
high main beam efficiency is desired, then much tighter tolerances 
are required. Even at X/50, which is usually considered to be a 
highly accurate antenna, the power loss due to errors is 6 percent. 

When the influence of errors on the side-lobe gain is considered, 
the results, as shown in Figure 10, depend strongly on the type of 
the error. If the error is an "overall" one (spherical aberration, 
cos 2c|), grating side lobes) , the relative side-lobe increment is 

• This quantity is the growth in side— lobe intensity expressed 
as a fraction of the degraded main-lobe intensity. Thus, even at 
a reflector-surface rms error of X/50. the growth in side-lobe peak 
power is only about 20 dB down from that of the main lobe. 
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Incidentally, the radiation pattern calculated for the "scallop" 
error shows the side-lobe gain to be decreased by the error. Where 
the power lost from the main lobe goes is a question - perhaps to 
the much higher lobes. At any rate, the geometrical error occurring 
in radial-rib types of antennas do not seem to present problems 
relative to the near-in side lobes. 

The "local" errors produce much lower effects if they are ran- 
dom. Thus, for example, if the "error-area" -to-aperture-area ratio 
-3 

IS 10 (a diameter ratio of 1/30) , then the reflector surface can 
have an rms error of X/16 and still have the side-lobe growth to be 
more than 30 dB down. 
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TABLE I. SECONDARY PATTERN CHARACTERISTICS PRODUCED 
BY A DISTRIBUTION (1 - r^)P OVER A CIRCULAR 
APERTURE WITHOUT PHASE ERRORS 


p 

G, GAIN FACTOR 

2 p + 1 
(p +1)^ 

0, HALF-POWER 
WIDTH 

POSITION OF 
FIRST ZERO 

FIRST SIDE LOBE, 
dB BELOW PEAK 
INTENSITY 

0 

1.00 

1.02 1 

. -1 1.22X 
sxn p 

17.6 

I 

0.75 

1.27 A 

. -1 1.63X 
sin 

24.6 

2 

0.56 

1.47 1 

. -1 2.03A 
sin ^ 

30.6 

3 

0.44 

1.65 1 

. -1 2.42X 
sin ^ 

- - - 

4 

0.36 

1.81 A 

1 

. -1 2.79X 
D 

- - - 
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Circular Aperture 



Figure 2. The effect of linear phase errors on the 

secondary pattern of a uniformly illuminated 
antenna . 




Figure 3. The effect of linear phase errors on the 
secondary pattern of an antenna with a 
quadratic illumination. 
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Figure 4. The effect of quadratic phase errors on the 

secondary pattern of an antenna with quadratic 
illumination . 


Figure 5. The effect of phase errors on the 

secondary pattern of a spherical reflector 






Figure 6. The effect of small-scale linear phase errors 
on the secondary pattern of a uniformly il- 
luminated antenna. 







151 








Relative 
side- lobe 
increment 


Square aperture, 
periodic checkerboard error 
(grating side lobe, = 45°) 


Ci rcular 
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Figure 10. Effect of errors on antenna side lobes. 
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APPENDIX A 

PROGRAM FOR SERIES EXPANSIONS 

The following computer program was used to compute the ex- 
pressions in Eq. (19) in the preceding text. In the present form, 
it computes F^. By changing Line 60 of the computer printout to 
4, it can be used to compute G^. In a similar manner, 5 computes 
H^, 6 computes I^, and 8 computes K^. The ranges of variables are 

n = 0, 1, 2, 3 , 4 

u - 1, ... 20 

For this program, use the following series of keystrokes; 


2nd 

X 

CMS 

STO 

06 

Stores 

value 

for U 

X 

STO 

11 

Stores 

value 

for n 

RST 

CLR 






R/S 


Use 


2nd 


FIX 


4 


the first time. 
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APPENDIX B 
PROGRAM FOR M AND N 


The following computer program expresses the integrands of 
the expressions for M(u) and N(u) in Eq. (41) of the preceding 
text. By using Label A, this program can be called up as part of 
a library program for numerical integration (using Simpson's rule). 
Changing the cosine term to a sine term in Step 014 relates the 
program to the integral for N(u) instead of M(u) . The digit 4 
in front of each integral has been included to simplify the en- 
suing computations . 
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SECTION 1 
INTRODUCTION 

The successful development of any structure depends largely 
on the identification of the critical or primary loads and other 
requirements on which the design is based. The future promises 
large structures which must be deployed, erected, assembled, or 
fabricated in space. For such structures, which will not be 
required to face the launch environment, the primary design re- 
quirements will be derived from the space-flight environment 
and will deal with phenomena as primary criteria which have been 
considered as only secondary in the past. The design of such 
genuine "space" structures will require a solid foundation of 
critical criteria. 

This chapter deals with a preliminary investigation of the 
tensions required to flatten metallized films or thin metal foils 
after they have been packaged in a folded condition. 

Metallized films and thin metallic foils are candidate 
materials for the surfaces of a variety of large space structures. 
The reflecting metal surface is usually required for solar or 
radio frequency reflectors where a smooth surface is desired for 
best performance. Normally, the material must be folded in the 
launch package and then deployed in space. It is of particular 
interest to determine the tension required to pull the material 
sufficiently flat for the mission. 

Astro Research Corporation has investigated previously the 
shortening effect of transverse folds in tapes (see ref. 1). A 
simplified model of a material with alternating folds in one 
direction is examined here. Also, experimental results on metal 
foil samples with both two- and three-dimensional folds are 
presented. 
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SECTION 2 

ANALYSIS OF ALTERNATING PARALLEL FOLDS 


Consider a length of material of thickness t with alternating 
parallel folds under a tension N as shown in Figure 1. Let s be 
the distance along the surface perpendicular to the folds and 2b 
the spacing between folds. Then for a unit width Z the moment on 
the film is 


M = NwJl 


and the curvature is given by 


d0 

ds 


iL 

El 


where 


cos 0 = 


dw 

ds 


Substituting Eqs. (1) and (3) into Eq. (2) gives 


cos 0 • d0 = wdw 

El 


For tensions where the surface is nearly flat 


( 1 ) 


( 2 ) 


(3) 


(4) 


0 = j atz=0,w=0 

Integrating Eq. (4) , we get 



0 


• d0 



(5) 
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which results in 


1 - sin 6 


2 

Nil w 
El 2 


( 6 ) 


or 


w 



sin 9 ) 


Calculating the root mean square displacement gives 



Substituting Eq, (6) and the following 


dz 


sin 9 
cos 9 


dw 


/W 

/I - sin 9 


into Eq. (7) results in 





(1 - sin 9) 


1/2 


sin 9 d(sin 9) 


(1 - sin'^ 9 ) 


1/2 


( 8 ) 


Letting y = sin 9 we get 




2EI 

Nil 


9 



1/2 


ydy 

/I + y 


(9) 
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Unfortunately, the model is not useful for a work-hardening mater- 
ial such as aluminum. 
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SECTION 3 

EXPERIMENTAL TESTS 


3.1 TESTS ON AN ALUMINUM FOIL STRIP 
WITH ALTERNATING PARALLEL FOLDS 


Recognizing that the analysis provided limited understanding 
of the tensions required to flatten a creased aluminum foil, three 
tests were conducted on creased aluminum foil samples. The specifi- 
cations of the test samples were: 


Aluminum . 

Thickness 

Width 

Length 

Length of folded section . 
Fold spacing 


Reynolds Wrap 
2.54 X 10"^ m 
0.10 m 
1.20 m 
1.00 m 
0.020 m 


broiling foil 


The samples were tested on a horizontal table dusted with 
talc to reduce friction. Tension was applied by dead weights via 
a string and pulley system. The length of the folded section, the 
height of the folds, and the half angle at the folds were measured 
under increasing load. A photograph of a test sample under load 
is shown in Figure 3. Results of these tests are shown in Figures 
4 and 5. The test data has been included in Appendix A. 

The fold angle 6 q at the tip of each fold proved to be very 
difficult to measure even with the load removed from the foil 
strip. And, because of the work hardening characteristic of the 
aluminum, this angle cannot be predicted by analysis. However, 

Eq. (6) has been plotted along with experimental data in Figure 5 
for an assumed constant 9 q of 60®. The experimental measurements 
of fold height shown on Figure 5 are somewhat ambiguous since at 
small loads the measurement represents 2 Wq, and at high loads the 
measurement is closer to w^. 
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3.2 TESTS ON A SQUARE FOIL SURFACE 

WITH FOLDS IN TWO DIRECTIONS 

A second series of tests was conducted on a square foil sample 
with the following dimensions; 

Thickness 2.54 x 10“^ m 

Sides 0.457 m 

Fold spacing ...... 0.020 m 

Arcs were cut in each edge of the foil, removing about 10 
percent of the area, and nylon cord catenaries were attached. 

The film was folded and packaged using the scheme shown in the 
photographs of Appendix B. Known loads were applied to the cate- 
naries in small increasing increments. The expansion of the film 
and average fold height were measured, and the radius of the cate- 
nary was estimated for each condition. Photographs of the expanding 
foil and test data are included in Appendix B. Results of the 
tests are shown in Figures 6 and 7. 
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SECTION 4 
CONCLUSIONS 

This preliminary investigation has shown that the effects of 
work hardening at creases in folded aluminvun foil require very high 
tensions in order to flatten the film in the region of the creases. 
Plastic effects also defeat the usefulness of analytical predic- 
tions. Fortunately, there are alternate surfaces available for 
both optical and radio frequency applications that large eliminate 
the difficulty associated with deploying foils. Optical surfaces 
can use polymer film coated with very thin layers of metal, mini- 
mizing the problems associated with foils and reducing the mass 
per unit area as well. Knit meshes are available for radio fre- 
quency reflector surfaces which, if they are tensioned, perform 
well and eliminate creasing in the packaged condition. 

Initially this preliminary study was viewed as a means of 
developing a model to predict the behavior of metallized films as 
well as metal foils. However, the difficulty experienced in pre- 
dicting the behavior of the "plastic hinge" at each crease, and 
in measuring the fold angles experimentally, eliminated any pos- 
sibility of extending the results. It appears that to establish 
the loads required to flatten metallized films, tests must be 
conducted on the specific coated film material. Plastic effects 
also occur when metallized film is creased, and although the prob- 
lems may be smaller than those experienced with foil, serious 
consideration should be given to packaging techniques which avoid 
folding or which minimize the number of folds. 
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Figure 1. Alternating parallel creases in a film or foil. 
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Figure 2. Membrane tension required to flatten 
creased film. 
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Figure 3 . 


Test sample under load. 


173 


Elongation under load. 



Edge tension, N/m 


Figure 4 . 


Extension of folded foil section under tension 



Fold height under load, 



Figure 5. Fold height of folded foil under tension. 
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Fold height, 


0.015 



0 20 40 60 80 

Deployment, % 


Figure 7. Deployment of a foil surface folded 
in two directions. 
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APPENDIX A 

TESTS ON ALUMINUM FOIL STRIPS WITH ALTERNATING PARALLEL FOLDS 

Three identical tests were performed on two separate occasions 
to study the effects folds have on deploying or unfolding a folded 
foil surface. The test consisted of a strip of aluminum foil 0.1-m 
wide and 1.2-m long. A series of 50 2-cm-wide folds were made in 
the center leaving 0.1 m at both ends for handling purposes (see 
Figures A-1 and A-2) . After all 50 folds were accurately made, the 
folded strip was placed in a vise and the folds set. The test was 
performed on a flat Formica table top dusted with talc powder to 
reduce friction. One end of the strip was firmly fixed by masking 
tape and the opposite end had a large paper clip anchor folded and 
applied so as to distribute the applied tension load evenly over 
the width of the strip (see Figures A-1 and A-3) . The loads were 
applied to a weight hanger on a string which ran over a pulley to 
the anchor on the aluminum strip. The length measured was from 
the first fold to the last fold. At each load case, the load was 
gently applied and the table top vibrated so as to allow the strip 
to overcome any friction in the system and extend to an equilibrium 
point. The length was measured and recorded as the loaded exten- 
sion (see Table A-1) . The load was then removed and the table 
vibrated to allow the strip to contract. The length, fold height, 
and fold angle were measured and recorded under load-relieved values 
(see Table A-1) . 

In the first test, starting at about Load Case No. 19 or 20 
and continuing on all subsequent load cases, an interesting phenom- 
enon occurred when the load was removed. The folds, instead of 
retracting in the manner it was originally folded, folded in the 
opposite direction (see Figures A-4 and A-5) . Unfortunately, it 
was not noticed until Load Case No. 21. This phenomenon probably 
accounts for the anomaly at the bottom of the percent relative 
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contraction versus percent elongation (loaded) graph shown in 
Figure A-6. After Load Case No. 24, a roller was used to flatten 
the foil to remove any remaining creases. This measured length 
was found to be approximately 0.0005 m shorter than the original 
value . 

The recorded data was used to calculate the load (in Newtons) , 
the percent extension (loaded), the percent relative contraction, 
percent extension (loaded) squared, and the effective load. These 
results are shown in Table A-2. Figures A-6 and A-7 provide a 
graphic presentation of selected data. 

Two additional tests were performed which were identical to 
the first 1-D loading aluminum foil fold test. The test was re- 
peated in order to obtain additional data not recorded during the 
first test and to check repeatability of the test. 

The test setup (see Figure A-3) and procedure were identical 
to the original test. A new strip of aluminum foil was prepared 
in order to avoid any possible strain-hardening effects of the 
first strip on the data. 

At each load case, the weight was gently lowered and the strip 
allowed to extend to an equilibrium point. The table was vibrated 
to allow the strip to overcome any friction. The loaded extension 
(first to last fold) , fold height, and angle were measured and re- 
corded (see Tables A-3 and A-4 and Figures A-4 , and A-8 through 
A-11) . The load was then removed and the strip allowed to contract. 
Again, the table was vibrated. The load-relieved extension, fold 
height, and angle were measured and recorded (see Tables A-5 and 
A-6) . 

At Load Case No. 8, it was found that the data needed on the 
fold angle was that of the fold tip and not that of the sides as 
was being recorded (see Figure 1 in main text) . From Load Case No. 

9 and on, the fold tip angle was recorded as noted with the data 
(see Tables A-4 and A-6) . 
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While lowering the weight at Load Case No. 15, the weight was 
accidentally dropped and pulled the strip further than the static 
load would have. Subsequently, the weight was increased to where 
the load extended the strip beyond the dynamic extension. As a 
result, the data at several loads previously planned for were not 
obtained. The test was continued as normal from then on. The 
recorded data was used to calculate the effective load, the percent 
extension (loaded) , the average fold angle (loaded) , and the average 
fold height (loaded) (see Table A-7) . As a result of these two 
abnormalities, it was decided to repeat the test. The same foil 
strip was refolded and tested again (see Table A-8) . 

In the second test, five preselected random folds were used to 
obtain data on the fold height and angle. It was originally hoped 
to obtain better representative data. Though the five folds did 
not necessarily represent an average of the 50 folds, some were 
very irregular and uneven. As a result, in the third test the 
fold height was the average of all the folds and the fold angle 
was obtained from that fold which best approaches the average fold 
angle (see Table A-8). In addition, the load-relieved fold height 
and fold angle were not recorded as was done in the second test 
because they were not needed. The recorded data was used to calcu- 
late the effective load and percent extension (loaded) (see Table 
A-9-) . Figure A-12 shows the fold height graphed against the per- 
cent extension at the last 10 percent of deployment using data 
obtained from the third test. 
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TABLE A-1. A1 strip, 1-D LOADING, FIRST TEST, RECORDED DATA: 

LOAD, LOADED EXTENSION, AND LOAD RELIEVED EXTENSION 
ANGLE AND HEIGHT 


LOAD 

CASE 

LOAD 

WEIGHT 

(g) 

LOADED 

EXTENSION 

(m) 

LOAD 

RELIEVED 

EXTENSION 

(m) 

LOAD RELIEVED 

ANGLE 

(“) 

HEIGHT 

(m) 

1 

0 

0.0095 

0.0095 



0.020 

2 

1.91 

0.1200 

0.0350 

— 

0.020 

3 

3.83 

0.1980 

0.0525 

4 

0.020 

4 

5.74 

0.2920 

0.0935 

6 

0.020 

5 

7.65 

0.3845 

0.1385 

7 

0.020 

6 

9.56 

0.4890 

0.2360 

15 

0.0195 

7 

11.48 

0.6585 

0.4185 

29 

0.0195 

8 

13.39 

0.7745 

0.5570 

36 

0.0150 

9 

15.30 

0.7925 

0.6175 

43 

0.0145 

10 

17.21 

0.8130 

0.6535 

48 

0.0140 

11 

19.13 

0.8300 

0.6780 

50 

0.0135 

12 

21.04 

0.8635 

0.7480 

57 

0.0115 

13 

22.95 

0.8695 

0.7540 

60 

0.0105 

14 

29.01 

0.8885 

0.7865 

63 

0.0095 

15 

42.72 

0.9170 

0.8465 

68 

0.0080 

16 

56.43 

0.9320 

0.8880 

72 

0.0070 

17 

70.14 

0.9410 

0.9095 

75 

0.0065 

18 

97.56 

0.9530 

0.9385 

85 

0.0055 

19 

152.40 

0.9660 

0.9585 

88 

0.0035 

20 

256.43 

0.9770 

0.9685 

88 

0.0025 

21 

407.24 

0.9845 

0.9735 

89 

0.0015 

22 

607.24 

0.9895 

0.9765 

* 

0.0010 

23 

1107.24 

0.9945 

0.9815 

* 

* 

24 

1542.24 

0.9965 

0.9860 

* 

* 


*Data not measurable for this case. 
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TABLE A-2. A1 STRIP, 1-D LOADING, FIRST TEST, CALCULATED FROM 
DATA: LOAD, PERCENT EXTENSION (LOADED), PERCENT 

RELATIVE CONTRACTION, PERCENT EXTENSION (LOADED) 
SQUARED, AND EFFECTIVE LOAD 


LOAD 

CASE 

LOAD 

(N) 

EXTENSION 

(LOADED)* 

(%) 

RELATIVE 

CONTRACTION 

(%) 

EFFECTIVE 

LOAD 

(N/m) 

EXTENSION 

(LOADED) 

SQUARED 

(%) 

1 

0.0187 

0.95 

0.0 

0.0 

0.009 

2 

0.0187 

12.00 

70.83 

0.187 

1.44 

3 

0.0376 

19.80 

73.48 

0.376 

3.92 

4 

0.0563 

29.20 

67.98 

0.565 

8.53 

5 

0.0750 

38.45 

63.98 

0.750 

14.78 

6 

0.0938 

48.90 

51.74 

0.940 

23.91 

7 

0.113 

65.85 

36.45 

1.130 

43.36 

8 

0.131 

77.45 

28.08 

1.310 

59.99 

9 

0.150 

79.25 

22.08 

1.50 

62.81 

10 

0.169 

81.30 

19.62 

1.69 

66.10 

11 

0.188 

83.00 

18.31 

1.88 

68.89 

12 

0.206 

86.35 

13.38 

2.06 

74.56 

13 

0.225 

86.95 

13.28 

2.25 

75.60 

14 

0.285 

88.85 

11.48 

2.85 

78.94 

15 

0.419 

91.70 

7.69 

4.19 

84.09 

16 

0.553 

93.20 

4.72 

5.55 

86.86 

17 

0.688 

94.10 

3.35 

6.90 

88.55 

18 

0.957 

95.30 

1.52 

9.55 

90.82 

19 

1.49 

96.60 

0.75 

14.9 

93.32 

20 

2.51 

97.70 

0.85 

25.1 

95.45 

21 

3.99 

98.45 

1.12 

39.9 

96.92 

22 

5.96 

98.95 

1.31 

59.6 

97.91 

23 

10.9 

99.45 

1.31 

109 

98.90 

24 

15.1 

99.65 

1.05 

151 

99.30 


*Lp = 0.9995 m (measurement taken after foil was rolled flat) 
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TABLE A- 3. Al STRIP, 1-D LOADING, SECOND TEST, RECORDED DATA 
LOAD, LOADED EXTENSION, AND FOLD HEIGHT (LOADED) 


LOAD 

CASE 

LOAD 

(8) 

LOADED 

EXTENSION 

(m) 

FOLD HEIGHT (LOADED) (m) 
(five random folds) 

FOLD #3 

FOLD #15 

FOLD #25 

FOLD #29 

FOLD #43 

1 

2.1 

■n 

0.0205 

0.020 

0.020 

0.020 

0.020 

2 

3.8 


0.021 

0.0195 

0.0195 

0.020 

0.0195 

3 

5.5 


0.020 

0.0195 

0.0195 

0.0195 

0.0195 

4 

10.5 


0.017 

0.0175 

0.017 

0.0175 

0.017 

5 

15.6 

0.735 

0.0125 

0.0130 

0.012 

0.0125 

0.0125 

6 

20.6 

0.846 

0.009 

0.010 

0.009 

0.009 

0.0095 

7 

25.6 

0.872 

0.0075 

0.0085 

0.0075 

0.008 

0.008 

8 

35.6 

0.907 

0.0055 

0.0060 

0.0055 

0.0055 

0.0055 

9 

49.6 

0.932 

0.005 

0.0045 

0.0040 

0.004 

0.004 

10 

63.8 

0.942 

0.004 

0.0035 

0.003 

0.0035 

0.003 

11 

77.9 

0.951 

0.0035 

0.002 

0.0025 

0.0025 

0.0025 

12 

92.0 

0.955 

0.0035 

0.002 

0.0020 

0.002 

0.0025 

13 

120.3 

0.961 

0.003 

0.002 

0.0015 

0.0015 

0.0025 

14 

162.9 

0.967 

0.0025 

0.002 

0.001 

0.001 

0.0015 

15 

653.2 

0.991 

0.0005 

0.001 

0.0005 

0.0005 

0.001 

16 

1524.6 

0.994 

* 

* 

* 

* 

* 

17 

2386.3 

0.998 

* 

* 

* 

* 

* 


*Too small to measure. 
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TABLE A-4. A1 STRIP, 1-D LOADING, SECOND TEST, RECORDED 
DATA: FOLD ANGLE (LOADED) 


LOAD 

n A CT7 

FOLD ANGLE (LOADED) (°) 
(five random folds) 








FOLD #3 

FOLD #15 

FOLD #25 

FOLD #29 

FOLD #43 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

15 

8 

10 

9 

0.0 

3 

16 

8 

10 

9 

8 

4 

40 

32 

34 

33 

35 

5 

60 

55 

60 

60 

60 

6 

72 

70 

75 

76 

73 

7 

80 

75 

78 

77 

75 

8 

85 

82 

85 

85 

82 

9 

23 

0 

37 

30 

10 

10 

25 

30 

30 

30 

15 

11 

25 

35 

30 

30 

20 

12 

25 

40 

32 

30 

20 

13 

28 

43 

35 

30 

20 

14 

28 

50 

45 

30 

25 

15 

60 

75 

75 

45 

68 

16 

65 

80 

90 

85 

80 

17 

90 

90 

90 

90 

90 


Note: Load Cases 1 through 8 show average angle; Load Cases 

9 through 17 show tip angle. Estimated accuracy ±5°. 
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TABLE A-5. Al STRIP, 1-D LOADING, SECOND TEST, RECORDED 

DATA: LOAD RELIEVED EXTENSION AND FOLD HEIGHT 

(LOAD RELIEVED) 


LOAD 

CASE 

LOAD 

RELIEVED 

EXTENSION 

(m) 

FOLD HEIGHT (LOAD RELIEVED) (m) 
(five random folds) 

FOLD #3 

FOLD #15 

FOLD #25 

FOLD #29 

FOLD #43 

1 

0.017 

0.020 

0.020 

0.020 

0.020 

0.020 

2 

0.027 

0.020 

0.020 

0.020 

0.020 

0.020 

3 

0.047 

0.020 

0.020 

0.0195 

0.0195 

0.0195 

4 

0.254 

0.0195 

0.019 

0.019 

0.019 

0.0195 

5 

0.521 

0.016 

0.0165 

0.016 

0.017 

0.017 

6 

0.711 

0.013 

0.0135 

0.0135 

0.0135 

0.013 

7 

0.760 

0.012 

0.0125 

0.0125 

0.013 

0.0125 

8 

0.854 

0.0115 

0.006 

0.010 

0.0105 

0.010 

9 

0.897 

0.010 

0.008 

0.0065 

0.009 

0.008 

10 

0.924 

0.009 

0.0025 

0.0055 

0.0055 

0.006 

11 

0.930 

0.0085 

0.0025 

0.005 

0.0065 

0.0055 

12 

0.938 

0.008 

0.0025 

0.004 

0.0055 

0.004 

13 

0.948 

0.0075 

0.0025 

0.003 

0.0045 

0.003 

14 

0.956 

0.0005 

0.001 

0.001 

0.0005 

0.0005 

15 

0.983 

* 

0.0005 

0.0005 

* 

* 

16 

0.983 

* 

0.0005 

* 

* 

* 

17 

0.987 

* 

* 

* 

* 

* 

18 

0.9995t 

* 

* 

* 

* 

* 


*Data not measurable for this case 
fMeasurement taken after foil was rolled flat. 
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TABLE A-6. A1 STRIP, 1-D LOADING, SECOND TEST, RECORDED 
DATA: FOLD ANGLE (LOAD RELIEVED) 




FOLD ANGLE (LOAD RELIEVED) (“) 


LOAD 


(five random folds) 


CASE 

FOLD #3 

FOLD #15 

FOLD #25 

FOLD #29 

FOLD #43 

1 

0.0 

0.0 

0.0 

0.0 

0.0 

2 

0.0 

4 

0.0 

5 

2 

3 

2 

4 

0.0 

6 

4 

4 

20 

14 

15 

15 

16 

5 

42 

40 

40 

38 

35 

6 

55 

55 

58 

58 

60 

7 

62 

55 

60 

62 

60 

8 

68 

73 

68 

65 

65 

9 

15 

5 

10 

30 

5 

10 

17 

40 

30 

30 

10 

11 

20 

40 

30 

30 

15 

12 

20 

40 

30 

30 

15 

13 

20 

42 

32 

30 

17 

14 

40 

60 

55 

50 

45 

15 

60 

68 

65 

50 

60 

16 

65 

70 

90 

50 

65 

17 

90 

90 

90 

90 

90 


Note: Load Cases 1 through 8 show average angle; Load Cases 

9 through 17 show tip angle. Estimated accuracy ±5°. 




TABLE A-7. Al STRIP, 1-D LOADING, SECOND TEST, CALCULATED 
FROM DATA: EFFECTIVE LOAD, PERCENT EXTENSION 

(LOADED) , AVERAGE FOLD ANGLE (LOADED) , AND 
AVERAGE FOLD HEIGHT (LOADED) . 


LOAD 

CASE 

EFFECTIVE 

LOAD 

(N/m) 

EXTENSION 

(LOADED) 

(%) 

AVERAGE 
FOLD ANGLE 
(LOADED) 

C) 

AVERAGE 
FOLD HEIGHT 
(LOADED) 
(m) 

1 

0.21 

4.2 

0.0 

0.020 

2 

0.37 

12.5 

8 

0.020 

3 

0.54 

20.8 

10 

0.0195 

4 

1.03 

50.4 

35 

0.017 

5 

1.53 

73.5 

60 

0.0125 

6 

2.02 

84.6 

73 

0.0095 

7 

2.51 

87.2 

77 

0.008 

8 

3.49 

90.7 

84 

0.0055 

9 

4.86 

93.2 

20 

0.0045 

10 

6.26 

94.2 

26 

0.0035 

11 

7.64 

95.1 

28 

0.0025 

12 

9.02 

95.5 

29 

0.0025 

13 

11.80 

96.1 

31 

0.002 

14 

15.98 

96.7 

53 

0.0015 

15 

64.06 

99.1 

65 

0.001 

16 

149.52 

99.4 

80 

<0.0005 

17 

234.02 

99.8 

<90 

<0.0005 


Note: 


Load Cases 1 through 8 show average angle; Load Cases 
9 through 17 show tip angle. Estimated accuracy ±5°. 
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TABLE A-8. Al STRIP, 1-D LOADING, THIRD TEST, RECORDED DATA: 
LOAD, LOADED EXTENSION, FOLD ANGLE (LOADED) , FOLD 
HEIGHT (LOADED), AND LOAD RELIEVED EXTENSION. 


LOAD 

CASE 

LOAD 

(g) 

LOADED 

EXTENSION 

(m) 

FOLD 

ANGLE 

(LOADED) 

(°) 

FOLD 

HEIGHT 

(LOADED) 

(m) 

LOAD 

RELIEVED 

EXTENSION 

(m) 

1 

10.5 

0.482 

4 

0.017 

0.263 

2 

20.5 

0.835 

7 

0.009 

0.664 

3 

30.5 

0.889 

8 

0.0065 

0.772 

4 

35.5 

0.903 

9 

0.0055 

0.803 

5 

40.4 

0.914 

9 

0.005 

0.827 

6 

45.3 

0.925 

12 

0.0045 

0.856 

7 

50.3 

0.930 

15 

0.004 

0.868 

8 

55.2 

0.935 

15 

0.004 

0.884 

9 

60.1 

0.939 

17 

0.0035 

0.893 

10 

70.1 

0.944 

20 

0.003 

0.907 

11 

80.1 

0.949 

22 

0.0025 

0.922 

12 

90.1 

0.953 

22 

0.002 

0.930 

13 

104.2 

0.957 

25 

0.002 

0.941 

14 

132.6 

0.962 

28 

0.0015 

0.950 

15 

175.2 

0.968 

30 

0.001 

0.960 

16 

237.0 

0.975 

35 

0.0005 

0.968 

17 

308.9 

0.980 

45 

0.0005 

0.973 

18 

422.6 

0.985 

50 

<0.0005 

0.977 

19 

536.2 

0.989 

60 

<0.0005 

0.979 

20 

763.4 

0.992 

60 

<0.0005 

0.981 

21 

1199.4 

0.995 

67 

<0.0005 

0.985 

22 

1628.1 

0.997 

67 

<0.0005 

0.989 

23 

2057.2 

0.9975 

80 

<0.0005 

0.992 

24 

2489.3 

0.998 

85 

<0.0005 

0.992 

25 

4978.6 

0.9985 

-~90 

<0.0005 

0.993 
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TABLE A-9. A1 STRIP, 1-D LOADING, THIRD TEST, CALCULATED 

FROM DATA: EFFECTIVE LOAD AND PERCENT EXTENSION 

(LOADED) 


LOAD 

CASE 

EFFECTIVE 

LOAD 

(N/m) 

EXTENSION 

(LOADED) 

(%) 

1 

1.03 

48.2 

2 

2.01 

83.5 

3 

2.99 

88.9 

4 

3,48 

90.3 

5 

3.96 

91.4 

6 

4.44 

92.5 

7 

4.93 

93.0 

8 

5.41 

93.5 

9 

5.89 

93.9 

10 

6.87 

94.4 

11 

7.86 

94.9 

12 

8.84 

95.3 

13 

10.22 

95.7 

14 

13.00 

96.2 

15 

17.18 

96.8 

16 

23.24 

97.5 

17 

30.29 

98.0 

18 

41.44 

98.5 

19 

52.59 

98.9 

20 

74.87 

99.2 

21 

117.63 

99.5 

22 

159.67 

99.7 

23 

201.75 

99.75 

24 

244.13 

99.8 

25 

488.25 

99.85 
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Figure A-4 . Load Case No. 14 (load relieved) 
(note reverse retraction) 

(Test No. 2) . 
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Normal 
retraction 
(load relieved) 



Figure A- 5. Illustration of reverse retraction which 
occurred at high loads when load relieved 
(see also Figure A-4) . 


193 







80 85 90 95 100 

Extension (loaded) , % 

Percent relative contraction vs. percent extension 
(loaded) (from data of first test) . 
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Figure A- 8. Load Case No. 4 (loaded) 
(Test No . 2 ) . 




Figure A-11. Load Case No. 15 (loaded) 
(Test No. 2) . 







Extension (loaded) , % 


Figure A-12. Fold height vs. percent extension 
(loaded) (from data of third test) 





APPENDIX B 

TEST ON A SQUARE FOIL SURFACE WITH FOLDS IN TWO DIRECTIONS 

The test sample used in these 2-D loading tests consisted of 
a square sheet of alximinum foil with 18.0-inch-long sides. The 
material is of the same type used in the 1-D loading tests (see 
Appendix A) . From this square, catenaries were cut in each side 
so the total area removed was 10 percent of the area of the original 
square (see Figure B-1) . A series of 0.25-inch-wide mylar tabs 
were fixed 2.0 cm apart on each of the four catenaries and a nylon 
cord was threaded through the tabs (see Figures B-2 and B-3) . A 
sleeve made of glass tape was wrapped around the two cords at each 
corner in order to keep the cords parallel to the catenaries at the 
corners . 

The manner in which the foil was folded was designed to best 
approach a typical folding pattern used in an actual application. 

The goal was to obtain a high-density, stowed package and a deploy- 
ment with no or little angular momentum being generated. This 
folding method will work equally well with a triangular surface. 

The folding guides were measured and creased by folding the square 
into a quarter section and then folding that over a straight edge. 
The folds were spaced 2.0 cm apart creating 10 folds between the 
center and the edge. Once the folding guides were made and the 
surface unfolded back into a square, the foil was folded as shown 
in Figures B-4 through B-7). 

The loads were applied to the catenaries by way of the cords. 
The final stowed configuration was placed on a table with a pulley 
at each corner over which the cords passed connected to weight 
hangers (see Figure B— 8) . For each load case, a small increment 
of weight was applied to the cords and the table was vibrated to 
remove any friction with the foil sheet or in the pulleys. After 
the foil had reached an equilibrium point, the length between 
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the center and each of the corners was measured Csee Table B-1) 
and a tracing of the outer perimeter of the surface was made. 

From this tracing, a planimeter was used to find the deployed area. 
From the data in Table B-1, the load, average length, and percent 
total length deployment were calculated (see Table B-2) . 

Upon initial loading, the "arms" unrolled to most of their 
full length before the main body of the foil surface unfolded (see 
Figures B-9 and B-10) . Beginning with Load Case No. 7, the foil 
surface had taken a shape so that an approximate radius of the 
sides could be obtained (see Figures B-11 through B-16) . Using 
the radius and the segment height, the area of the segment could 
be calculated. From this, the deployed surface area is found and 
averaged with the area using the planimeter (see Table B-3) . 

After the final load case. No. 17, a roller was used to flatten 
the foil to remove any remaining creases and a tracing was made. 

A second test was performed which was similar to the first 
2-D loading test to obtain information on fold/wrinkle height not 
recorded during the first test. The test setup was similar to that 
used in the first test except larger pulleys with ball bearings 
were used and the test was performed on the surface table (see 
Figure B-17) . The original sheet of aluminum foil was refurbished 
by replacing about 25 torn tabs and three out of the four corners. 
Before folding, the foil sheet was rolled flat, but still contained 
some 0.005-inch-high wrinkles from its previous use. The sheet 
showed signs of wear and as a result it was not possible to load 
it to the higher loads used in the first test. 

At each load case, the weights were gently applied and the 
distance from each corner to the center measured and recorded (see 
Table B-4) . To measure the fold/wrinkle height, a platform mea- 
suring 5.45 inches by 4.45 inches was set at a height such that 
when placed under the sheet it slightly pushed it up. This height 
was accurately measured. The platform was placed between two 
corners with one edge along the side of the sheet. A clear plexi- 
glass disc was placed on the foil sheet above the platform. The 
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height was measured at two points 1 inch to either side of center 
and 2.75 inches in from the side. The thickness of the disc was 
subtracted and the height recorded (see Table B-5) . This was 
performed at two opposite quadrants (see Figure B-18) . 

While folding the sheet, it was decided not to coil up the 
four "arms" at the last stage of folding (see Figures B-6 and B-7) . 
This was done in order to reduce wear and tear. As a result of 
this, while deploying the surface, an opposing fold pattern developed 
perpendicular to each side. The coiling would have allowed the foil 
surface to deploy without creating the opposing fold pattern. The 
opposing folds were removed by hand and the test was continued. 

It was originally intended to use the center— to— corner length 
to correlate the percent of deployment. Though the length shows 
little correlation to the load when comparing it to the first test, 
it was decided to use the load to calculate the percent of deployment 
(see Table B— 6) . A possible explanation of this abnormality with 
the length is that the corners are not the natural corners of the 
aluminum foil sheet. The corners are made with mystic Kapton tape. 
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TABLE B-1. SQUARE FOIL SURFACE, 2-D LOADING, RECORDED 
DATA: TENSION LOAD AND CENTER-TO-CORNER 

LENGTH 


I 













TABLE B-2. SQUARE FOIL SURFACE, 2-D LOADING, CALCULATED 

DATA: LOAD, AVERAGE LENGTH, AND PERCENT TOTAL 

LENGTH DEPLOYMENT 


LOAD 

LOAD, 2T 

AVERAGE 

LENGTH, 

DEPLOYMENT 
LENGTH 
(% total) 

CASE 

GRAMS 

NEWTONS 

Lavg 

(in. ) 

1 

0.0 

0.0 

2.25 

17.67 

2 

5.5 

0.054 

2.34 

18.38 

3 

19.8 

0.194 

2.56 

20.11 

4 

48.4 

0.475 

3.55 

27.89 

5 

77.0 

0.755 

5.44 

42.73 

6 

132.2 

1.297 

8.00 

62.84 

7 

245.6 

2.409 

10.05 

78.95 

8 

359.0 

3.521 

10.93 

85.86 

9 

472.4 

4.633 

11.39 

89.47 

10 

699.2 

6.857 

11.64 

91.44 

11 

926.0 

9.081 

11.71 

91.99 

12 

1358.3 

13.321 

11.85 

93.09 

13 

2265 

22.21 

11.99 

94.19 

14 

4600 

45.11 

12.13 

95.29 

15 

6875 

67.4 

12.24 

96.15 

16 

9150 

89.7 

12.30 

96.62 

17 

18200 

178.5 

12.51 

12.54* 

12.73** 

98.27 


*Measurement taken after foil was rolled flat. 
**Original length before test. 
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TABLE B-3. SQUARE FOIL SURFACE, 2-D LOADING, CALCULATED 
FROM DATA; EFFECTIVE LOAD, AREA, RADIUS, AND 
PERCENT TOTAL AREA DEPLOYMENT 


LOAD 

CASE 

EFFECTIVE 

LOAD, 

N = T/R 
(N/m) 

AREA 

(in2) 

RADIUS 

(in.) 

AREA 

DEPLOY^IENT 
(% total) 

1 

n/a 

n/a 

n/a 

n/a 

2 

n/a 

2.39 

n/a 

0.82 

3 

n/a 

2.69 

n/a 

0.93 

4 

n/a 

2.26 

n/a 

0.79 

5 

n/a 

3.01 

n/a 

1.03 

6 

n/a 

14.82 

n/a 

5.08 

7 

5.50 

54.95 

8.61 

18.83 

8 

6.51 

93.2 

10.65 

31.96 

9 

6.44 

142.3 

14.17 

48.80 

10 

5.27 

209.1 

25.6 

71.71 

11 

6.00 

222.5 

29.8 

76.30 

12 

6.64 

240.4 

39.5 

82.44 

13 

9.02 

255.7 

48.5 

87.69 

14 

10.23 

264.7 

54.7 

90.78 

15 

23.52 

270.5 

56.5 

92.76 

16 

32.05 

274.2 

55.1 

94.03 

17 

63.08 

281.0 

55.7 

96.36 



287.4* 

58.91* 

98.56* 



291.6** 

60.19** 

100.0** 


^Measurement taken after foil was rolled flat. 
**Original length before test. 
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TABLE B-4. SQUARE FOIL SURFACE, 2-D LOADING, FOLD /WRINKLE 
HEIGHT TEST, RECORDED DATA; TENSION LOAD AND 
CENTER-TO-CORNER LENGTH. 


LOAD 

CASE 

TENSION 
LOAD, 2T 
(g) 

CENTER-TO-CORNER LENGTH 

(in.) 





1 

680.4 

12.15 

12.45 

12.25 

12.00 

2 

1247.4 

12.25 

12.55 

12.45 

12.20 

3 

2827 

12.35 

12.70 

12.45 

12.25 

4 

4536 

12.40 

12.60 

12.50 

12.60 

5 

6804 

12.35 

12.75 

12.50 

12.75 
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TABLE B-5. SQUARE FOIL SURFACE, 2-D LOADING, FOLD/WRINKLE 

HEIGHT TEST, RECORDED DATA: FOLD/NRINKLE HEIGHT 



FOLD /WRINKLE HEIGHT (in.) 

LOAD 

CASE 

QUADRANT 

BETWEEN 1 AND 2 

QUADRANT BETWEEN 3 AND 4 



h2 

h3 


1 

0.433 

0.561 

0.446 

0.511 

2 

0.376 

0.366 

0.266 

0.408 

3 

0.178 

0.131 

0.236 

0.274 

4 

0.138 

0.111 

0.134 

. 0.133 

5 

0.100 

0.080 

0.081 

0.074 
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TABLE B-6. SQUARE FOIL SURFACE, 2-D LOADING, FOLD/WRINKLE 

HEIGHT TEST, CALCULATED FROM DATA: TENSION LOAD 

AVERAGE LENGTH, AVERAGE FOLD/WRINKLE HEIGHT, 
PERCENT TOTAL AREA DEPLOYMENT BASED ON LOAD AND 
EFFECTIVE LOAD. 
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Foil surface being folded 

(note 2.0-cm-wide folding guides) 






















Figure B-17 . Fold/wrinkle height test setup 
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CHAPTER 4 


METEOROID DAMAGE TO ROD-TYPE STRUCTURAL ELEMENTS 


by 

John M. Hedgepeth 
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INTRODUCTION 


Large space structures are likely to be built up from many 
rod-type structural elements. While the change of meteoroid damage 
to any particular element is very small, the chance of damage some- 
where in the structure may be appreciable because of the large 
total length of the elements. Little is known of the meteoroid 
damage characteristics of linear elements besides the early treat- 
ment in Reference 1. No experimental data exist at present. This 
report contains further theoretical study of the problem. 

METEOROID STATISTICS 

Let n (m) be the density function per unit volume for meteoroids 
of mass m and let the velocity components of each meteoroid be ran- 
dom variables with Gaussian statistics. Let the probability dis- 
tribution function be denoted as 


f (V) 


n (m) 




(v^+v^+v^)/2U^ 
\ X y zr 


( 1 ) 


where U is the rms velocity magnitude. Note that we have assiimed 
stationarity and isotropy. 

The expected number of meteoroids in a volume S with a mass 

between m and m+dm and velocity components in the intervals 

(v^,v +dv ) , (v ,v +dv ) , (v ,v +dv ) is 
X X X y y y ' ' z z z 


dq = Sf (V) dm dv dv dv 
^ ' ' X y z 


which can also be written as 

= Sf(v) dm sin 0 dVd0d<t> 


dq 


( 2 ) 



where we have set 


V = V sin 9 cos 9 

X 

V = V sin 9 sin 9 

y 

V = V cos 9 
z 


Note 


2 ^ 2 ^ 2 

X y z 


= V‘ 


Flux 

Consider the area A in the x-y plane. We would like to 
determine the number of meteoroids of mass between m and m+dm with 
flow through the area A from one side per unit time. 



Select the lower side as the one through which the meteoroids are 
passing. Then all the meteoroids in the cylindrical volume shown 
above that have spherical-coordinate velocity components of V, 9, 
and (j) pass through the area A in the unit time. Thus the total 


223 



rate of flow from all directions coming from below is 


« j[/2 2 IT 


dm 


■// / 

0 0 *^0 


AV cos 9 f (V) sin 0 d0d<{idV 


( 3 ) 


Substituting from Eq. (1) and performing the integrations yields 


dQ _ n (m) U 


dm 




(4) 


Let N (m) be the number density of meteoroids of mass m or greater 
so that 


oo 

N (m) = J" n 

m 


(m) dm 


Let (m) be the number flux per unit area from one side of meteor- 
oids with mass greater than or equal to m. Then 


(m) 


UN(m) 

/ 6 ¥ 


(5) 


In many cases, the data for meteoroid frequency are approxi- 
mated by the formula 


N(m) = Km ^ 


( 6 ) 


Hits on a Long Cylinder 

Consider a circular cylinder of diameter d and length £ >> d 
as shown in the following figure; 
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z 



In the same way as before, all the meteoroids with velocity- 
component coordinates of V, 9, (j) in the above volume of V£d sin 0 
will hit the cylinder in a unit time. All meteoroids with velocity- 
component coordinates of V, 6, <j) outside of that volume will miss 
the cylinder. Thus, the total rate of hits from all directions is 

2tt 

V£d sin 9 f (v) sin 0 v^d(})d0dv 



Substituting for f (v) and integrating gives 


dQ 

dm 


TTdJl 


n (m) y 


(7) 
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The number of hits per unit time with a mass equal to or 
than m is 


H(m) 


00 


I 


m 


dm 


dm 


or 

H (m) = TTdJl $ (m) 


where we have made use of the flux equation, Eq. (5) . 


Average Meteoroid Velocity 

The "average" velocity ^ can be obtained as 

ave 


.«> ^ir TT 


''ave • Wh J f f ® 


0 0 -TT 


47T 




O 

/ 


^3 ^-3V^/2U^ 


Evaluating the integral gives 


V 

ave 



FRACTURES BY THE "SWEPT- AREA" CRITERION 

According to Reference 2, the penetration depth is, 
infinite targets 


P = k m°-352 pl/6 v2/3 

CO CO 


greater 


(8) 


(9) 


for semi- 


( 10 ) 
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where 


m = meteoroid mass, g 

= meteoroid density, g/cc 

V = meteoroid velocity, km/s 

= penetration depth, cm 

= a constant (= 0.42 for aluminum) 

According to Reference 3, Denardo, et al, find that the crater 
diameter is approximately twice the penetration. Let us assume 
that the long cylinder will fracture whenever the cylinder of 
radius equal to (Eq. (12) ) about the meteoroid path encompasses 
the entire cylinder. Thus, if any part of the cylinder is not 
swept out by the destructive zone, the cylinder is assumed to be 
unfractured. This criterion seems to be particularly appropriate 
for unidirectional filamentary composites. 

If P<„(ni,V) < d/2, then there is no fracture. If P^(m,V) > d/2, 
then the thickness of the volume containing meteoroids that will 
produce fracture is 2P^ - d. 



Therefore, the number of fractures per unit time will be (see 
Eq. (6)) 
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,oo ^7T _7r 


Nf - / am II /_ (2P^-d)V£ sin 9 f (v) sin 6 V^d<})d0dv 


0 V* •'0 -IT 


OO 00 


I 

0 


= TT^ I dm I (2P -d)V^J,f(V)dV 


( 11 ) 


where V* is the velocity for which P^ = d/2. Reversing the order 
of integration and substituting for f(V) yields 


N. 


2U' 


a 


(2P^-d)V^ g-3V^/2U^ n(in)dmdV 


(12) 


0 •'m* 


where, now m* is the mass (as a function of V) for which P^ = d/2. 
Thus 


.1/0.352 


m* (V) = 

\2k 


\ «> m 

From Eq. 

(6) we get 

n (m) = 

YKm~^“^ 


(13) 


Then integrating gives 


N, 


3 . . 33, ,2k pl/6 


2U' 



3 y 2 y-0.352 


C 

I 


^3+(2y/3 0.352) -3V^/2U^ 


dV 


(14) 
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Substituting t for 3V^/(2U^) 


gives 



(y/0.352)-! J 


^(y/1.056)+1 



( 15 ) 


where is the simple-minded estimate obtained by assuming that 

all meteoroids with larger mass than m* (U) (Eq. (13)) will fracture 
the cylinder if they strike it. Thus from Eqs. (5) and (8) 


N 


f 


0 



dHUN(m*[u]) 


(16) 


The integral in Eq. (15) can be evaluated in terms of the 
gamma function to yield 



( 2 / 3 )^/^*°^^ 

(y/0.352)-1 


■(- 


056 



(17) 


APPLICATION TO EARTH ORBIT 


For earth orbit, the meteoroid flux is given in Reference 4 to 


be 


$ = 4.26 X 10 nm 


ioT the mass m in grams greater than 10 The factor ri accounts 

for the effects of earth shadowing and defocusing. It varies from 
0.5 at the Earth's surface to a maximum of 0.691 at an altitude of 
three Earth radii and decreases to about 0.619 at GEO. We there- 
fore set, with only slight conservatism. 


n = 0.7 
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The average velocity is given as 20 km/s. Thus, from Eq. (9) 


U = 21.7 km/s 

Referring to Eq. (5), we can solve for N (m) to give 
N(m) = 5.97x10“^^ particles/m^ 

Thus, 

K = 5.97x10 particles/m^ 

Y = 1.213 

For this value of y 


N 


f _ (2/3) 


1.149 


N. 


2.45 


r (3.149) 


= 0.59 


The value of m* (U) is obtained from Eq. (13) to be 


m* (U) = 


[ 


2 X 0.42 X (0.5)^/^ X (21.7)^/^ 


] 


1/0.352 


where we have used k^ = 0.42 (aluminum) for want of a better value, 
and p = 0.5 g/cc . 

This gives 


m*(U) = 0.0067 d 


2.84 


where d is in centimeters. 
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Substituting into Eqs. (15) and (16) with = 1 m gives the 
fracture rate/m/s of 

= 0.59 X X ^ X 21,700 x 5.97 x lO”^® 

X (0.0067)"^*^^^ X (j“2. 84 X 1.213 

= 2 . 4 X 10 d ^ fractures/m/s 

= 7.6x10 fractures/m/yr 

For members 1 cm in diameter and a total length of 1 km, the proba- 
bility of a fracture in 10 years is about 0.008. On the other 
hand, if the diameter is only 1 mm and the length is 1 km, then 
the probability of a fracture in one year is about 0.2, which is 
too high. 


REFERENCES 


1. MacNeal, R.H.; Meteoroid Damage to Filamentary Structures. 

NASA CR-869, September 1967. 

2. Meteoroid Damage Assessment. NASA SP-8042, May 1970. 

3. Denardo, B.P.; et al. : Projectile Size Effects on Hypervelocity 

Impact Craters in Aluminum. NASA TN D-4067, 1967. 

4. Meteoroid Environmental Model - 1969 (Near Earth to Lunar 

Surface). NASA SP-8013, March 1969. 


231 



1. Report No. 

NASA CR-165744 

4. Title and Subtitle 


2. Government Accession No. 


Considerations In the Design of Large Space Structures 


3. Recipient's Catalog No. 


5. Report Date 
August 1981 


6. Performing Organization Code 


7. Author(s) 

John M. Hedgepeth, Richard H. MacNeal, Karl Knapp, and 
^ Ch arles S. MacGillivray 

9. Performing Organization Name and Address 

Astro Research Corporation 
6390 Cindy Lane 
Carpinteria, CA 93013 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


8. Performing Organization Report No. 
ARC-TN-1101 


10. Work Unit No. 

11. Contract or Grant No. 

NASl-15347 

13. Type of Report and Period Covered 
Contractor Report 

14. Sponsoring Agency Code 


I 15. Supplementary Notes 

NASA Langley Technical Monitor: Dr. Melvin S. Anderson 


16. Abstract 


This report contains several analytical studies of topics relevant to the 
design of large space structures. Topics covered are: 

1. The types and quantitative evaluation of the disturbances to 
which large Earth-oriented microwave reflectors would be subjected 
and the resulting attitude errors of such spacecraft 

2. The influence of errors in the structural geometry of the 
performance of radiofrequency antennas 

3. The effect of creasing on the flatness of tensioned reflector 
membrane surfaces 

4. An analysis of the statistics of damage to truss- type structures due 
to meteoroids . 


17. Key Words (Suggested by Author(s) | 

Large space structures 


18. Distribution Statement 


Unclassified - Unlimited 


19. Security Classif. (of this report) 
Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21. No. of Pages 22. Price 
231 


N-305 


For sale by the National Technical Information Service, Springfield, Virginia 22161 



